NASA/TP-2009-214643 



Further Development of Ko Displacement Theory 
for Deformed Shape Predictions of Nonuniform 
Aerospace Structures 


William L. Ko and Van Tran Fleischer 
NASA Dryden Flight Research Center 
Edwards, California 



September 2009 


NASA STI Program ... in Profile 


Since its founding, NASA has been dedicated 
to the advancement of aeronautics and space 
science. The NASA scientific and technical 
information (STI) program plays a key part in 
helping NASA maintain this important role. 

The NASA STI program operates under the 
auspices of the Agency Chief Information 
Officer. It collects, organizes, provides for 
archiving, and disseminates NASA’s STI. The 
NASA STI program provides access to the NASA 
Aeronautics and Space Database and its public 
interface, the NASA Technical Report Server, 
thus providing one of the largest collections of 
aeronautical and space science STI in the world. 
Results are published in both non-NASA channels 
and by NASA in the NASA STI Report Series, 
which includes the following report types: 

• TECHNICAL PUBLICATION. Reports of 
completed research or a major significant 
phase of research that present the results of 
NASA Programs and include extensive data 
or theoretical analysis. Includes compilations 
of significant scientific and technical data and 
information deemed to be of continuing 
reference value. NASA counterpart of peer- 
reviewed formal professional papers but has 
less stringent limitations on manuscript length 
and extent of graphic presentations. 

• TECHNICAL MEMORANDUM. Scientific 
and technical findings that are preliminary or 
of specialized interest, e.g., quick release 
reports, working papers, and bibliographies 
that contain minimal annotation. Does not 
contain extensive analysis. 

• CONTRACTOR REPORT. Scientific and 
technical findings by NASA-sponsored 
contractors and grantees. 


• CONFERENCE PUBLICATION. Collected 
papers from scientific and technical 
conferences, symposia, seminars, or other 
meetings sponsored or co-sponsored 

by NASA. 

• SPECIAL PUBLICATION. Scientific, 
technical, or historical information from 
NASA programs, projects, and missions, 
often concerned with subjects having 
substantial public interest. 

• TECHNICAL TRANSLATION. English- 
language translations of foreign scientific 
and technical material pertinent to 
NASA’s mission. 

Specialized services also include creating 
custom thesauri, building customized databases, 
and organizing and publishing research results. 

For more information about the NASA STI 
program, see the following: 

• Access the NASA STI program home page 
at http://www.sti.nasa.sov 

• E-mail your question via the Internet to 
help@sti.nasa.gov 

• Fax your question to the NASA STI Help 
Desk at 443-757-5803 

• Phone the NASA STI Help Desk at 
443-757-5802 

• Write to: 

NASA STI Help Desk 

NASA Center for AeroSpace Information 

7115 Standard Drive 

Hanover, MD 21076-1320 



NASA/TP-2009-214643 



Further Development of Ko Displacement Theory 
for Deformed Shape Predictions of Nonuniform 
Aerospace Structures 


William L. Ko and Van Tran Fleischer 
NASA Dryden Flight Research Center 
Edwards, California 


National Aeronautics and 
Space Administration 

Dryden Flight Research Center 
Edwards, California 93523-0273 


September 2009 


The structural shape prediction method described in this report is protected under U.S. Patent No. 7,520, 176, 
issued April 21, 2009. Therefore, those interested in using the method should contact the NAS A Innovative 
Partnership Program Office at the Dryden Flight Research Center for more information. 


Cover art: NASA Dryden Flight Research Center, art number 090014. 


NOTICE 

Use of trade names or names of manufacturers in this document does not constitute an official endorsement 
of such products or manufacturers, either expressed or implied, by the National Aeronautics and Space 
Administration. 


Available from: 

NASA Center for AeroSpace Information 
7115 Standard Drive 
Hanover, MD 21076-1320 
(443) 757-5802 


TABLE OF CONTENTS 


ABSTRACT 1 

NOMENCLATURE 1 

INTRODUCTION 2 

BASICS OF THE KO DISPLACEMENT THEORY 3 

NONUNIFORM CANTILEVER BEAMS 4 

Discretization 4 

Slope Equations 5 

Deflection Equations 6 

FURTHER MATHEMATICAL DEVELOPMENT 6 

Slope Equation 7 

Deflection Equations 7 

Limit-Case Mathematical Problems 8 

SLIGHTLY NONUNIFORM CANTILEVER BEAMS 9 

Slope Equations 9 

Deflection Equations 10 

UNIFORM CANTILEVER BEAMS 12 

Slope Equations 12 

Deflection Equations 12 

Relative Magnitude of Strain Terms 13 

TWO-POINT SUPPORTED BEAMS 15 

Left End Simply Supported 15 

Left End Fixed 17 

TWO-LINE STRAIN-SENSING SYSTEM 17 

Sensor Installation Examples 17 

Wing Box 18 

Cantilever Tubular Beam 18 

Deflections 18 

Cross-Sectional Twists 19 

APPLICATION TO FREE-FREE BEAMS 19 

APPLICATION TO STRUCTURAL PANELS 20 

VALIDATION OF KO DISPLACEMENT THEORY 20 

Input Strains 20 

Tapered Tubular Beams 21 

Cantilever Beams 21 



Bending Only 21 

Bending and Torsion 22 

Two-Point Supported Beams 23 

Simply Supported Ends 23 

Fixed Ends 24 

Tapered Wing Boxes 24 

Depth-Tapered Wing Boxes 25 

Unswept Wing Boxes 25 

Swept Wing Boxes 26 

Width-Tapered Wing Boxes 27 

Double-Tapered Wing Boxes 27 

EXACT AND EXPANDED DEFLECTION EQUATIONS 28 

CONCLUDING REMARKS 29 

FIGURES 31 


APPENDIX A: DERIVATION OF DEFLECTION EQUATIONS FOR NONUNIFORM CANTILEVER BEAMS... 58 
APPENDIX B: EXPANSION OF LOGARITHMIC TERMS 64 


APPENDIX C: DERIVATIONS OF DEFLECTION EQUATION FOR SLIGHTLY NONUNIFORM 
CANTILEVER BEAMS 67 

APPENDIX D: DERIVATIONS OF DEFLECTION EQUATIONS FOR UNIFORM CANTILEVER BEAMS 74 

APPENDIX E: RELATIVE MAGNITUDES OF STRAIN TERMS FOR UNIFORM CANTILEVER BEAMS 77 

REFERENCES 79 


IV 



ABSTRACT 


The Ko displacement theory previously formulated for deformed shape predictions of nonuniform beam 
structures is further developed mathematically. The further-developed displacement equations are expressed 
explicitly in terms of geometrical parameters of the beam and bending strains at equally spaced strain- 
sensing stations along the multiplexed fiber-optic sensor line installed on the bottom surface of the beam. 
The bending strain data can then be input into the displacement equations for calculations of local slopes, 
deflections, and cross-sectional twist angles for generating the overall deformed shapes of the nonuniform 
beam. The further-developed displacement theory can also be applied to the deformed shape predictions of 
nonuniform two-point supported beams, nonuniform panels, nonuniform aircraft wings and fuselages, and 
so forth. The high degree of accuracy of the further-developed displacement theory for nonuniform beams is 
validated by finite-element analysis of various nonuniform beam structures. Such structures include tapered 
tubular beams, depth-tapered unswept and swept wing boxes, width-tapered wing boxes, and double-tapered 
wing boxes, all under combined bending and torsional loads. The Ko displacement theory, combined with 
the fiber-optic strain-sensing system, provide a powerful tool for in-flight deformed shape monitoring of 
unmanned aerospace vehicles by ground-based pilots to maintain safe flights. 

NOMENCLATURE 

c half depth of uniform beam (distance from neutral surface to bottom or top surface), in. 

CG center of gravity 

t'i half depth of nonuniform beam at strain-sensing station, X( (distance from neutral surface 

to strain-sensing station, X [ ), in. 

c n value of q at free end (beam tip) strain-sensing station, x n , in. 

c n / co beam tip-root depth ratio 

c{x) half depth of nonuniform beam at point x, in. 

co value of q at fixed end (beam root) strain-sensing station, xq , in. 

d chord-wise distance between two span-wise parallel strain-sensing lines, in. 

deg degree 

dj chord-wise distance between tapered front and rear strain-sensing stations at { X( , x\ }, in. 

E Y oung modulus, lb/in 2 

e calculated nonzero deflection error at right fixed end of two-point supported beam, in. 

G shear modulus, lb/in 2 

/ moment of inertia, in 4 

i = 0, 1, 2, 3, ..., n, strain-sensing station identification number 

j index 

l length of beam, in. 

M bending moment, in-lb 

Mi bending moment at strain-sensing station, X ( , in-lb 

n index associated with the last span-wise strain-sensing station 

P applied load, lb 

SPAR Structural Performance And Resizing 

T applied torque, in-lb 

t thickness of wall, in. 

w width of uniform wing box, in. 

w n width of wing box tip, in. 

wq width of wing box root, in. 



x, y Cartesian coordinates (x in beam axial direction, y in vertical direction), in. 

X[ axial coordinate (or symbol) associated with i -th strain-sensing station, in. 

>7 deflection in y-direction at X [ , in. 

D 

yf deflection in y-direction at x/ for two-point supported beam, in. 

y(x) deflection in y-direction at axial location, x, in. 

a half of angle, sustained by d, separating two strain-sensing lines on tubular beam, deg 

/3 wing box swept angle, deg 

Yi surface shear strain, rad 

A / = l/n , distance between two adjacent strain-sensing stations, in. 

£j surface bending strain at strain-sensing station, X [ , in/in. 

£(x) surface bending strain at axial location, x, in/in. 

9j slope of deformed beam at strain-sensing station, Xj , rad or deg 

0 (x) slope of deformed beam at axial location, x, rad or deg 

V Poisson’s ratio 

( 7 • axial stress at strain-sensing station, x; , lb/in 2 

<7(x) surface bending stress at axial location, x, lb/in 2 

< pi cross-sectional twist angle at axial location, x; , rad or deg 

( )' quantity associated with rear strain-sensing line 

INTRODUCTION 

Conventional aerospace instrumentation systems, such as those required for in-flight deformed shape 
monitoring of aircraft, are subject to strict weight and size restrictions. The existing electro-optical flight 
deflection measurement system (ref. 1 ) and conventional strain gage measurement system are too heavy 
(because of strain gage lead wires) for lightweight flying vehicles. The most attractive candidate for an in- 
flight deformed shape measurement system is the fiber Bragg grating (FBG) sensing system, because the 
sensing lines are extremely small and lightweight (approximately the size of human hairs), possess excellent 
flexibility, and can be highly multiplexed with hundreds of sensors located on single fibers at desired 
sensing intervals. 

When these highly multiplexed sensors are installed on beam-like structures (such as aircraft wings), the 
measured strain data can be input into special algorithms to calculate the deflections and cross-sectional 
twist angles of the beam-like structure for deformed shape predictions. Special algorithms based on the Ko 
displacement theory are described in detail in reference 2. In the mathematical formulations, the beam is 
first discretized into multiple domains. Based on the classical beam theory, the displacement theory (ref. 2) 
was formulated to predict the deformed shapes (slopes, deflections, and cross-sectional twist angles) of 
cantilever beams (uniform, tapered, slightly tapered, stepwise tapered) using the fiber-optic measured strain 
data as inputs. 

The Ko displacement equations, combined with the onboard fiber-optic strain-sensing system (called 
Method For Real-Time Structure Shape-Sensing, U.S. Patent No. 7,520,176, issued April 21, 2009; ref. 3), 
provide a powerful tool for in-flight deformed shape monitoring of unmanned aerospace vehicles by 
ground-based pilots to maintain safe flights. In addition, the wing shape that is monitored in real time could 
then be input into the aircraft control system for aeroelastic wing shape control. 
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In the original formulation of the Ko displacement theory (refs. 2, 3) for nonuniform cantilever beams 
(for example, tapered, slightly tapered), the displacement equations are written in indicial recursion 
relationships because of multiple beam domains, and they need to be fully developed into more compact 
mathematical forms. In this report, the displacement equations for nonuniform beams (cantilever and two- 
point supported beams) formulated previously (refs. 2, 3) are further developed into more compact 
mathematical forms so that the beam deflections at each strain-sensing station can be expressed exclusively 
in terms of inboard geometrical parameters of the beam and inboard surface bending strains. The accuracy 
of the further-developed displacement theory for nonunifonn beams is validated by finite-element analysis 
of various nonuniform beam structures (tapered tubular beams, such as cantilever and two-point supported 
beams, depth-tapered unswept and swept wing boxes, width-tapered wing boxes, and double-tapered wing 
boxes) under different loading conditions. 

BASICS OF THE KO DISPLACEMENT THEORY 

The formulation of the displacement theories is built upon the classical beam differential equation 
(elastic curvature of a deformed uniform beam) given by (refs. 4, 5) 


d 2 y _ M{x) 
dx 2 El 

in which y is the vertical deflection, x is the span-wise coordinate, M(x) is the bending moment, E is the 
Y oung modulus, and / is the moment of inertia. 

At the cross section of the beam, x, the bending moment, M(x) (eq. (1)), induces the bending stress, 
<7(x) , at the bottom (or top) fiber of the beam (outermost point from the neutral axis) (refs. 2-5) with the 
magnitude proportional to M(x) according to 


a(x) = 


M(x)c 

I 


( 2 ) 


in which c is the half depth (distance from the neutral axis to the bottom or top surface) of the uniform 
beam. The bending stress, o(x) , is related to the bending strain, e(x ) , at the same point on the bottom (or 
top) surface of the beam through Hooke’s law as 


£(x) = 


<7(x) 

E 


( 3 ) 


When equations (2) and (3) are combined, the bending moment, M(x), can be expressed in terms of the 
bending strain, £(x) , as 


M(x) = El 


£(x) 

c 


(4) 
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In light of equation (4), beam differential equation (1 ) becomes 


= (5) 
dx z c 

Note that under the present strain formulation, beam differential equation (5) contains only the half depth of 
the beam, c , and the bending strain, £(x) . The flexural rigidity, El, is eliminated. 

Beam differential equation (5) for the uniform beam could be used with sufficient accuracy for 
nonuniform beams for which the cross section is assumed to change gradually (ref. 4, p. 143). Namely, the 
half depth of the beam, c, is no longer constant but is a function of x (that is, c = c(x) ). Hence, for the 
nonuniform beam, equations (4) and (5) can be modified, respectively, to 


M(x) = El 


£(x) 

c(x) 


( 6 ) 


and 


d^y = £(x) 
dx 2 c(x) 

The formulation of the Ko displacement theory (ref. 2) for nonuniform beams is based upon modified 
beam differential equation (7). By discretizing the beam into multiple domains, and assuming { c(x),e(x) } to 
be piecewise linear, one could integrate equation (7) over each domain to yield local slope and deflection 
equations for generating the overall deformed shape of the beam (ref. 2). The Ko displacement theory 
requires further mathematical development as described in the subsequent sections. 

NONUNIFORM CANTILEVER BEAMS 

The nonuniform cantilever beam is defined as a cantilever beam in which the cross section varies slowly 
along the axis of the beam. This type of beam can also be called a weak nonuniform beam. 

Discretization 

Figure 1 shows a nonuniform cantilever beam with a length, /, and varying depth, 2c(x) (for a tubular 
beam, c(x) is the local radius). Let the fiber-optic strain-sensing line be multiplexed with n +1 
(0, 1, 2, 3, ..., n) number of equally spaced strain-sensing stations, Xj (i = 0, 1, 2, 3, ..., n ), at axial 
locations, x = Xj , on the bottom surface of the beam, and let the measured bending strains be denoted by £; 
(i = 0, 1, 2, 3, ..., n). Note from figure 1 that the first and last strain-sensing stations, { xq,x„ }, are located 

at the fixed end (x = 0) and at the free end (x = Z), respectively. The index, /' (= 0, 1, 2, 3, ..., n), is also 
used as the identification number of the strain-sensing stations. The nonuniform cantilever beam is now 
divided into n sections with identical domain length, Al = l / n . The nonuniform beam can then be idealized 
as a piecewise linearly tapered beam (either tapering down or tapering up) between any two adjacent 
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strain-sensing stations, { Xj_ j , x,- }. Under this assumption, the half depth of the beam, c(x ) , in the region, 
X;_i <x<Xj, between the two adjacent strain-sensing stations, { X(_\ , x ( - }, can be expressed as a linear 
function of x as 


c{x) = c i _i-{c i _i-c i ) 1 1 ; Xj_ ] < x < Xj (8) 

Al 

in which { c,_] , c ; - } are the values of c(x) at the strain-sensing stations, { x,-_ | , x/ }, respectively. Likewise, 
the bending strain, £(x), in the region, x ( -_] < x < x/ , between the two adjacent strain-sensing stations, 
{ x;_] ,x ; - }, can be assumed to change linearly with x as 


£(x) — 1 1 


u) 


X-Xj-l 

Al 


x i-\ <x<Xj 


( 9 ) 


in which { £,_i,£ ; - } are the values of £(x) at the strain-sensing stations, { x/_ | , x ( - }, respectively. 

If c(x) is a weak linear function of x (small slope term in eq. (8)), then the bending moment for the 
nonuniform beam, M(x) = Ele(x)/ c(x) (eq. (6)), will vary almost linearly between the two adjacent strain- 
sensing stations, { x,-_ | ,x ; - }, because the higher order terms in the binomial expansion of the denominator, 

[c(x)] _1 , can be neglected. This argument may not hold, however, if c(x) is not a weak linear function ofx, 
because the higher order terms in the binomial expansion of the denominator [c(x)j 1 must be retained. 

For the uniform beam case, [ c(x) = c = constant], the bending moment, M(x), is directly proportional to 
the bending strain, £(x) , according to equation (4). Therefore, if the strain, £(x) , is a linear function of x, 
then M(x) is also a linear function ofx. 

Figure 2 shows the slopes, 9j (i = 0, 1, 2, 3, ..., n), and deflections, y/ (i = 0, 1, 2, 3, ..., «), at the 
strain-sensing stations, x,- (i = 0, 1, 2, 3, ..., n), on the lower surface of the deformed tubular cantilever 
beam (using the uniform tubular beam as an example). 

Slope Equations 

The slope, tan0(x) , of the nonuniform beam at the axial location, x, between the two adjacent strain- 
sensing stations, x/_| < x < x,- , can be obtained by integrating equation (7) once and enforcing the 
continuity of slope at the next inboard strain-sensing station, x ( -_ | . Namely, 


jc 

tan 6(x) = [ - — ^-dx 

Jx i - 1 dx 2 

Slope increment 


+ tan 9j_i 
Slope at x- t _ | 


f> + 

J U K(x) * 


X/-1 < X < Xi 


( 10 ) 


in which tan0/_i is the slope at x\_\ . 
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Substituting beam half depth equation (8) and strain equation (9) into slope equation (10), and carrying 
out the integration to the upper limit, x = x; , causes equation (10) to yield the slope, tan 0, [ = tan Q(x l ) J, at 
the strain-sensing station, Xj , as (see Appendix D of ref. 2) 


tan 9 1 = Al 


£ i — 1 £ i . £ i—\ c i £ i c i — 1 . c i 

c i—l~ c i (C;_j — Cj) c i— 1 


+ tan (9;_| 


( 11 ) 


When i = l, the slope, tan 9\_\ = tan 9q , will be the slope at the left fixed end, which is zero ( tan 9q = 0 ) for 
the cantilever beam. If the reference left end is simply supported (as in the two-point supported beam case), 
then the slope, tan0q > is nonzero ( tan (9q ^ 0 ). 


Deflection Equations 

The deflection, y(x) , of the nonuniform beam in the region, x,_] <x<X[, between the two adjacent 
strain-sensing stations, {xj_i,Xj}, can be obtained by integrating slope equation (10) and enforcing the 
continuity of deflection at the adjacent inboard strain-sensing station, x,_j . Namely, 


rx 

y(x)= tan 9{x)dx + y,--l 

J xi_i — * — ' 

- — - — - — — - — - Deflection 

Integration of slope a t x •_ \ 


£ ( x ) , , f* „ 

ax ax + tan t/,-_ i ax + y, _i 

,x i-l Jx i-1 c(x) J N-I * n ^4. 

Deflection 


(X rx 

J X: i J X: 


Deflection increment Deflection at x 
due to tand ( _] 


at Xj_i 


( 12 ) 


(x,-i < X < Xj ) 


in which y;_i is the deflection at X;_j . 

Carrying out the double and single integrations of deflection equation (12) in light of beam depth 
equation (8) and strain equation (9), we obtain the deflection, y,- [ = y(x, ) ], at the strain-sensing station, x, , 
as (see Appendix D of ref. 2) 


)’i = (MY 


£ i - 1 £ i £ i-\ c i £ i c i — 1 


2( c i-[-c i ) ( C; _, -c ; )- 


Ci 


c i^Og e — — + (Cj-l-C/) 
c i - 1 


• + >’,■_! + Altm 6 j_\ 


(13) 


When i = l, the deflection and slope, { y\_\ = yq , tan0j_i = tan0q }, become the deflection and slope at the 
reference left fixed end where ' yq = 0. tan 0q =0 }. If the left end is simply supported (as in the two-point 
supported beam case), the deflection, yq , is zero (yq = 0 ), but the slope, tanOq , is nonzero ( tan0q z 0 ). 

FURTHER MATHEMATICAL DEVELOPMENT 

Equations (11) and (13) are the recursion formulae. In this section these equations are further developed 
to more compact mathematical forms expressed explicitly in terms of beam geometry, {Al , q } , and 
bending strains, £ ; - (fig. 1). 
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Slope Equation 


Equation (11) is a recursion formula showing that the slope, tan 9j , at the strain-sensing station, x h is 
expressed in terms of the slope, tan (9;_i , at the adjacent inboard strain-sensing station, .17 _ | . Applying the 
descending indices relationship causes slope equation (11) to become the summation form as 
(see Appendix A for details) 


tan 6 1 = Al ^ 

7= 1 


£ j-i £ .i ! £ j~\ c j £ j c j~ 1 1 () „ c j 

c j- 1 ~ c j (c j- 1 - Cj ) 2 c ;-l 


+ tan 9q 


(14) 


in which tanOq is the slope at the left fixed end (reference strain-sensing station, xq) and is zero 
( tan (9 q ^0 ). For a simply supported left end (as in the two-point supported beam case), however, the slope, 
tanOq, is nonzero ( tan 9q ^0). Note that the slope, tan 9j , is determined from the geometrical and strain 
data evaluated at all the inboard strain-sensing stations (xq , xj , X2 , . . . , x ; - ), including the current strain- 
sensing station, x, . 


Deflection Equations 

Under the present discretization formulation, deflection equation (13) is also a recursion formula 
showing that the deflection, y ( - , at the strain-sensing station, x ; - , is expressed in terms of the deflection, 
>7_l , and the slope, tan 0/_i , at the adjacent inboard strain-sensing station, x ( _| . By applying the 
successive descending recursion relationships for >’,■_] in light of equation (13), and writing equation (14) 
for tan 0/_i , we could express deflection equation (13) for y t with two series summations (with different 
summation limits) containing only the geometrical parameters and bending strains (see Appendix A for 
derivations). Namely, 


7=1 


£ j - 1 £ i £ j-i c j £ j c j - 1 


2 (c ; -l-c ; ) (Cj-l-Cj) 3 


Cj log* 


C M 


+ (Cj_l~Cj) 


Contributions from deflection terms 


i — 1 




7=1 


£ J-1 £ j , £ j~ ]C j £ j c j~ 1 


- + 


c j'-l c i 


J 


(Cj-l-Cj) 


loge 


J ' 


c i - 1 


(15) 


+ V() + (i)Al tan 9q 


Contributions from slope terms 


in which yq = tan 9q = 0 at the fixed end (reference strain-sensing station, xq ). For a simply supported left 
end (as in the two-point supported beam case), however, the slope, tan (9 q , is nonzero ( tan (9 q ^0). Note 
that the deflection, y,- , is expressed in terms of the geometrical and strain data evaluated at all the inboard 
strain-sensing stations ( X(y , x\ , X2 , . . . , x ; - ), including the current strain-sensing station, X[ . 
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After grouping the terms, we can rewrite equation (15) in an alternative form (see Appendix A for 
details) as 


>7 = (M 2 £ 

7= 1 

i— 1 


[2(i -;) + l] 


£ /-i £ j £ ;'-i c ; £ j c j - 1 


2 (c ; -_! c ; ) ( C j_ { - C j ) 3 


Cj l°g e 


_^7_ 

C J -1 


+ (Cj_l~Cj) 


+ {Al) 2 J j {i-j) | 
7=1 


£ 7-1 £ 7- £ 7 £ 7— l lQge ^/_ 


(Cy-i-Cy) 


7' 


c /-l 


(16) 


+ V() + (i)Al tan 6q 


in which yp = tan 6q = 0 at the fixed end (reference strain-sensing station, xq ). 

By inputting the geometrical parameters of the beam and the bending strain data into deflection equation 
(16), we can calculate the deflection of the beam, y,- , at any strain-sensing station, jq . Therefore, the overall 
deformed shape of the nonuniform cantilever beam can be generated without the need to know the flexural 
rigidity of the beam, EL The ability to generate deformed shapes without this information is a powerful 
characteristic of the Ko displacement theory, because it does not require the tedious and time-consuming 
calculations of the moment of inertia 7(x), which changes along the span of a complex beam-like structure 
such as an aircraft wing. 

Note that the fixed-end condition, yp = tan Op =0 , will cause the calculated deflection, >7 (or slope, 
tan 0j ), to be relative to the tangent line passing through the reference strain-sensing station, xq , at the fixed 
end. Thus, deflection equation (16) is applicable to both stationaiy and moving cantilever beams such as 
aircraft wings during flight. 

Furthermore, deflection equation (16) is formulated for the nonuniform cantilever beam under lateral 
loading only (bending). If the beam is also subjected to axial loading, the measured strains, £[ 

(i = 0, 1, 2, 3, ..., n), will contain both bending and axial-load strain components. To remove the axial-load 
strain components, an additional strain-sensing line is needed on the upper surface along the axis of the 
beam. If the neutral axis is located at the half depth of the beam, the axial-load strain components can be 
removed by averaging the magnitudes of the lower and upper surface strains measured at the same strain- 
sensing cross section to yield the true magnitude of the bending strains. 

Limit-Case Mathematical Problems 

The mathematical forms of slope equation (14) and deflection equation (16) cannot be applied directly 
to the uniform beam, because at the limit of (q_j = c/), the logarithmic terms and the denominators with 
(c;_i-c/) factors in equations (14) and (16) will go to zero (that is, log e (cq/<q_i ) = 0 , q_j-c/= 0), 
causing equations (14) and (16) to be indeterminate. The way to circumvent this limit-case mathematical 
problem is to first consider the case of a slightly nonuniform cantilever beam and use the perturbation 
method to expand the logarithmic terms in series form. The resulting series expansion will contain 
(c/_| -<q) factors in the numerators, which will cancel out the (q_j-c/) factors in the denominators. 
Thus, the modified mathematical forms for equations (14) and (16) will be able to degenerate into the limit- 
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case expressions for uniform beams, ( q_] = q). The modified slope and deflection equations for slightly 
nonuniform cantilever beams are derived in the subsequent sections. 

SLIGHTLY NONUNIFORM CANTILEVER BEAMS 

When the cantilever beam is slightly tapered (that is, (q/q_i) — » 1), the logarithmic terms in the slope 
and deflection equations (eqs. (14) and (16)) will approach zero (that is, log e (c//q_j)— >0). Therefore, the 
perturbation method must be used to expand the logarithmic term, log e (q/q_i) , in the neighborhood of 
(c//c ; -_| ) ~ I to obtain nonzero mathematical expressions, so that the slope and deflection equations can 
degenerate into those for the uniform cantilever beam case, (c ; -/q_| ) = 1 . 

Slope Equations 

For a slightly tapered beam, the logarithmic term, log e (q/q_i) , in slope equation (11) can be 
expanded in the neighborhood of (q/q_| ) = 1 (that is, log e (q /q _| ) ~ 0 ). Carrying out the expansion up to 

9 

the second-order terms in (c,-_ j -q) causes the logarithmic term, log e (q7q_i) , to take on the following 
form as c, approaches q_ ] [(q/q_i) — » 1 ] (see Appendix B for details of mathematical expansions; 
refs. 2, 6): 


log e 


c i 
c i - 1 


2cU 


c i 


c i - 1 


*1 


(17) 


Substitution of equation (17) into slope equation (11) yields the following slope equation for the slightly 
nonuniform cantilever beam (ref. 2), with the ( Cj_\ - c/ ) factor in the denominators eliminated: 


4/ 

7 x 

C; 





tan 0; = 

2q-t. 

2 1 
c i-\ / 

£/-!+£/ 

+ tan 0,_| 

; q_, c f 

(18) 


Applying the descending recursion relationship causes slope equation (18) to become 


n Al i 1 
tan 0i= — £ l 


2 ~ Cj - 1 


f 


c j 


2 -- 
v c j- [ 


£ j~l +£ j 


+ tan 0 O ; q-i^q 


(19) 


Note that slope equations (18) and (19) do not contain the ( q_] - Cj ) factor and can now be used for the 
limit case (q_ j = q) of the uniform beam. 
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Deflection Equations 

When the beam becomes slightly tapered (that is, q — » q_i), the term, [q log e (q/q_i ) + (q_| -q)] > 
in equation (13) can be expanded in the neighborhood of (q/q_i ) ~ 1 (that is, log e (q/q_| ) ~ 0 ). Carrying 
out the expansion up to the third-order terms in (q_] -q) causes the term, 
[q log e (q/q_i) + (q_j -q)] , to take on the following form as q approaches c l -_ 1 [(q/q_i)— > 1 ] (see 
Appendix B for details of mathematical expansions; refs. 2, 6): 


qlog e ' +(q_ i q) 

(q — l -q ) 2 

r 3+ (c i _i-c,)'i 

L c 7-l J 

6q_i 

c i - 1 


c i 


< 7- 1 


*1 


In light of equation (20), deflection equation (13) becomes 


>7 


(Air 


6 c, 


7—1 


Vf \ 


3 C «- 

A-\ +£ i 

_V c i- 1 j 



+ >7-1 + Al tan0,_| 


Ci 


Ci - 1 


*1 


( 20 ) 


( 21 ) 


which contains no ( q_ j - q ) factor in the denominators. 


Applying the descending recursion relationship of equation (21) for y/_| , and writing slope 
equation (19) for tan0/_i causes deflection equation (21) to take on a form that has two series summations 
(see Appendix C), as 


>7 = 


(AD 


2 i 


1 




3-- 


V 


c 7-l 


£,•_!+£/• 


| (Al) 2 ^ (i-j) 


J 


(=1 C J ~ 1 



f 

X 

- 


2- 

c j 

£ j-l +£ j 


c j~ l J 


V 



Contributions from deflection terms 


Contributions from slope terms 


+ yq + (i)Al tan 0q 

( 22 ) 


in which yq = tan 6q = 0 at the fixed end (reference strain-sensing station, xq ) of the cantilever beam. Note 
from equation (22) that the two summation upper limits are different. 

Through grouping similar terms in the two series summations, we can write equation (22) in a final 
compact form consisting of only one series summation (see the derivation in Appendix C). Namely, 


37 



3(2; -1)- (3j-2)- 


c i~j + 1 
c i~j 


£i- j + (3 j - 2 )£i_j+i \ + >'() + (i)Al tan0 o (23) 


In equation (23), vq = tan 0q = 0 for the cantilever beam; therefore, the terms, yq + (i)Al tan 0q , can be 
dropped. The yq + (i)Al tan (9q terms are purposely kept, however, to make equation (23) applicable to other 
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beam cases, such as simply supported beams, for which {yp = 0. tan 9q ^0} at the reference strain-sensing 
station, xq (see the section entitled, “Two-Point Supported Beams”). Note that equation (23) is applicable 
to both stationary and moving cantilever beams (for example, airborne aircraft wings), because the 
deflection, y ( - , calculated from equation (23) is relative to the tangent line ( tan 0p = 0) stemming from the 
reference strain-sensing station, xq , which can be either stationaiy or moving (fig. 2). 

For easy visualization of the functional progression with ascending indices, equation (23) can be written 
out explicitly for different indices, /' ( = 1,2,3, ..., n), as 


fl = 


(Al) z 


6c 0 


IY 3 

3- — 

£ 0 +£l 

A c o) 



+ y>Q+ Al tan 0q 


(23a) 


T2 =- 


(AlY 


f>[ 

Y 3 

3 — — 

£l +£ 2 

1 

+ — 

Y 3 

9-4 — 

£p+4£l ] 


A q ) 

. 

c 0 

A c 0 ) 

JJ 


+ V() + 2 Al tan 9q 


(23b) 


>’3 = - 


{Aiy 


1 

y 7 

3- — 

£ 2 + £3 

1 

+ — 

Y 3 

9-4 — 

£\ + 4 £ 2 

1 

+ — 

Y 3 

15-7 — 

£q+4£i 

1 

A 2 

A C 2J 


ci 

A c i ) 


CO 

V cpj 


J 


+ yp + 3 Al tan 0p 


(23c) 


F4 


f>[ 

Y 3 

~ 

1 

+ — 

Y 3 

- 

1 

+ — 

Y 3 

~ 

3-^ 

£ 3 +£4 

9—4 — 

£ 2 +4 £ 3 

15-7 — 

£\ + 2 £ 2 

l C 3 

A c 3 ) 


C 2 

A c 2 ) 


Cl 

V c \ ; 



+ 


+ - 


c 0 


21-10 — 
c 0j 


£q -t-lOq 


■ + yp + 4A1 tan 9q 


(23d) 


T5 


1 

Y 3 

- 

1 

+ — 

Y 3 


1 

+ — 

Y 3 


3-^ 

£4+ £5 

9-4 — 

£ 3 + 4£ 4 

15-7 — 

£ 2 + 7£ 3 

c 4 

.V C4 y 


C3 

A c 3 ) 


c 2 

V c 2) 



1 

+ — 

q 


Y 3 

21-10 — 

+ 10^2 

1 

+ — 

Y 3 

27-13 — 

£q+13£i ] 

A c i ) 


CO 

A c oJ 

JJ 


+ yp + 5AI tan 0p 


(23 e) 


y n 


(Al) 


2 n 


" 1 
I — 


6 /=! c n-j 


3(2;-!)- (3; -2) 


c n-j 


£ n -j + (3 j - 2)£ n - j+ i \ + >’p + nAl tan Qq (23f) 


in which the fixed-end condition, yp = tan0p = 0 , holds. 

Equation (23) explicitly shows that the deflection, y,- , at the strain-sensing station, X [ , is the function of 
the inboard half depth of the beam ( cp , q , C 2 , q ) and the associated inboard strains ( £p , q , £2 
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including the values of { q , £/ } at the current strain-sensing station, Xj , in which the deflection, > 7 , is 
evaluated. By inputting the geometry of the beam and the measured strain data into deflection equation (23), 
we can calculate the deflection, >7 , at any strain-sensing station and thereby generate the deformed shape of 
the slightly nonuniform cantilever beam such as an aircraft wing. Note that deflection equation (23), 
developed for the slightly nonuniform beam, can now be used for the uniform beam case by setting 
q = q_ ] = c as a limit case without creating the mathematical breakdown problem encountered in 
deflection equations (15) and (16). 

UNIFORM CANTILEVER BEAMS 

When the nonuniformity of the beam diminishes (that is, q = q-_ | = c ; fig. 2), slope equations (18) and 
(19) and deflection equations (21) and (23), developed for the slightly nonuniform cantilever beam, can 
degenerate into the limit forms for the uniform cantilever beam. 

Slope Equations 

For the limit case of a uniform beam, q = q_ j = c, slope equations (18) and (19) developed for the 
slightly nonuniform cantilever beam degenerate, respectively, into the following limit forms for the uniform 
cantilever beam: 


tan 6/ = 


AI_ 

2c 


(£,■_ 1 +£,■) + tan 0 f _! 


(24) 


Al 1 

tan 6j = — X i£ j-\ + £ j 1 + tan e () 

j = 1 


(25) 


Deflection Equations 

Also, as the nonuniformity of the beam diminishes (that is, q = q_ | = c), deflection equations (21) and 
(23) for the slightly nonuniform cantilever beam degenerate, respectively, into limit forms for the uniform 
cantilever beam as 


yi = 


{Al ) 2 

6c 


(2 Ei-\ + £j ) + > 7-1 + Al tan 0;_i 


(26) 


and 


/ A 7 l 

yi = X [( 3 -/ - ] { £ i-j + { 3 J - 2 ) £ i-j+\ ] + >'()+ (i)Al tan 0 o 
7=1 

in which yq = tan = 0 holds at the left fixed end. 


( 27 ) 
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After grouping the strain terms of identical indices, we can modify equation (27) into the following 
more compact form (see details of derivation in Appendix D): 


yi = 


(Al) 2 

6 c 


i — 1 

(3/ - l)£q + 6 ( i - j)sj + E[ 

7=1 


+ y f ) + (i)Al tan Qq 


(28) 


Equation (28) can be written out explicitly for different indices, i (= 1, 2, 3, n), so that the functional 

progression of equation (28) with the increasing value of i can easily be seen: 


vi 


(Alf 

6c 


(2eq + ) + yq + Al tan 6q 


(28a) 


yi 


(Al) 2 

6c 


(5£q + 6e\ + +£2 ) yq + 2 Al tan 0q 


(28b) 


T3 = 


(Al) 2 

6c 


[8£q + 6(2e\ + £2 ) + £3 ] + yq + 3 Al tan 0 q 


(28c) 


>’4 


(Al) 2 

6c 


|^1 1£q + 6(3£j + 2 £2 + £3 ) + £4 J + yq + 4 Al tan 0 q 


(28d) 


>’5 


(Al) 2 

6c 


[14£ 0 + 6 ( 4 £ ] + 3£2 + 2£3 + £4 ) + £3 + yq + 5 Al tan Oq 


( 28 e) 


y n = 


(Al) 2 

6 c 


n — 1 

(3 n - 1)£ 0 +6 X ~ 7') £y + £ n 

7=1 


+ yq + 11 Al tan 0() 


(28f) 


If the terms, yq and Al tan 6q , are removed, then deflection equation (28f) will be identical to deflection 
equation (12e) of reference 2 for the uniform cantilever beam. Equation (28) explicitly shows that the 
deflection, y,- , at the strain-sensing station, jq , is obtained by summing up only the inboard strains 
(£q ,£i, E j , ■■■,£{ ), including the strain, £; , at the current strain-sensing station, Xj . 

Relative Magnitude of Strain Terms 

The tip deflection, y n , expressed in equation (28f) is the summation of the strain data obtained at all the 
inboard strain-sensing stations. This section examines the relative magnitude of those strain terms, because 
each strain term has a different coefficient. The example used is the classical uniform cantilever beam under 
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tip vertical load, P. For this case, the bending strain, £(x), is a linearly decreasing function of the axial 
coordinate, x, and becomes zero at the beam tip (refs. 4, 5). Namely, 


£(x) = 


( 

1 -- 
V 1 J 


£0 


(29) 


in which l is the length of the beam, and £q i s the maximum bending strain at the fixed end ( x = xq ). Based 
on the strain function given in equation (29), the beam tip deflections, y n (eq. (28f)), for different values of 
n can be calculated as follows (see Appendix E for details): 

For n = 1 : 


(A/) 2 , l 2 

y \ - — — (2^0 + £ \ ) - ~ £ 0 

6c 3c 


(30a) 


For n = 2: 


(Aiy x 

yi - — — (5£o + 6£i + £ 2 ) 

be 


l_ 

3c 


£0 


f 

CF625 +0.375 + 0 
V maximum j 


A l 2 


3c 


£0 


(30b) 


For n = 4: 


y 4 - ^ - - [11 £q +6(3^ + 2£ 2 +£3) + £4] 


6c 

2 ( 


l_ 

3c 


£0 


0.34375 + 0.421875 + 0.1875 + 0.046875 + 0 


^ / 2 


V 


maximum 


J 


3c 


£0 


(30c) 


For n = ; 


>’8 


(A IT 

6c 

,2 f 


[23£q +6^7 £j + 6£ 2 + 5£3 + 4£4 + 3£^ + 2£g + £7 ) + £g J 




A 


0.17969+ 0.28711 +0.21094 + 0.14648 + 0.09375 + +0.05273 + 0.02344 + 0.00586 + 0 
V maximum j 


3c 


£0 


(30d) 


In equations (30b)-(30d), the dominant terms are underlined and labeled, “maximum.' 
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Figure 3 compares the relative magnitude of the strain terms, inside the brackets of equations 
(30a)-(30d), for different values of n. For n = 1-2, the maximum strain term is located at the fixed end ( xq ). 
When n is greater than 2 (n >2 ), the maximum strain term shifts toward the next outboard strain-sensing 
station, x \ , which becomes the most dominant strain-sensing station. Thus, the accuracy of measured strain 
at the dominant strain-sensing station greatly affects the accuracy of the calculated beam deflections using 
equation (28). 


TWO-POINT SUPPORTED BEAMS 

For the case of a nonuniform cantilever beam, the fixed end (yq = tan 6q = 0 ) is located at the left-hand 
side at the reference strain-sensing station, xq . For a two-point supported beam (simple beam), the two 
ends can be under any support condition (fixed, simply supported, or elastically supported). 

Left End Simply Supported 

For a simple beam with a simply supported left end, the deflection is zero (yo = 0), but the slope is 
nonzero (tandq^O). If deflection equation (23) (yq = tan 0q=O) for the cantilever beam is used to 
calculate the deflection curve of the simple beam case, the deflection curve will have a zero slope 
( tan (9 q = 0 ) at the left support, and a maximum deflection ( y n Z 0 ) at the right end (deflection curve at 
position A, fig. 4). Thus, to apply cantilever beam deflection equation (23) for the simple beam case, the 
unknown slope, tandq - at the left support point must be determined first. By setting the right-end deflection, 
y n , to zero ( y n = 0 ) in equation (23f), we can determine the unknown slope, tan (9q , at the simply 
supported left end as 


_ Al v 1 
tanflq =- — X — 


6/7 Cj-{ 


3(2j — 1) — (3y — 2) - 


c n—j + 1 
c n-j 


£ n—j (3./ y+l [ — n/ \y n (31) 


in which the y n is the deflection of the cantilever beam tip given by 


y n = 


(Al) 


2 n 


n 1 

I- 


6 Z'J - > 


3(2;-l)-(3j-2)- 


c n—j + 1 
c n-j 


£ n- j (3./ 2)s n —j + 1 


(32) 


Equation (23-f); vq =tan 9q =0 


which is equation (23f) with the terms, yq +nAl tan dq , removed, because yq = tan (9q = 0 for the cantilever 
beam. 
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With tan0Q determined in equation (31), deflection equation (23) for the cantilever beam can be 

D 

modified to calculate the deflection, y, , of the simple beam with a simply supported left end (the right end 
can be under any support condition) as 


B W 

yi =-r- 


j= 1 C J-1 


3(2 j -1)- (37-2)- 


Ci- 


i~j + 1 
c i~j 


£i—j+(3j 2 


Yj [equation (23); _V()=tan6()=0] 


i 

■ -y n 

n 

Shift factor 


(33) 


In equation (33), the last term, (i / n)y n ( i = 1, 2, 3, ..., n), is a linearly increasing shifting factor (with 
increasing /') for vertical downward shifting of the cantilever beam deflection curve from position A to 
position B, to convert it into the simple beam deflection curve (fig. 4; zero shift at the left end). Namely, the 
downward shifting is equivalent to the clockwise rotation of the inclined reference line OA, with respect to 
the left support point 0, by 9q deg to the horizontal position OB with the length slightly shrunk (fig. 4). For 
easy visualization of the functional progression with increasing i, equation (33) is written out for different 
strain-sensing stations as 


b (Air i 
yi = —^\ — 

6 q 


B (Mf 1 
F3 =^^) — 
6 c 2 


VA 

V c i) 


y i 


£ 1 +£ 2 


b m~ 


6c 0 

J_ 

c 0 


( \ 
3- — 

£ 0 +£l 

l c 0 J 



-y n 


+ - 


V 4 £f 


V 


c oJ 


£ 0 +4£i 


-y n 


( \ 
3 — — 

£ 2 +£3 

1 

+ — 

Y v 

9-4 — 

£j + 4£ 2 

1 

+ 

Y v 

15-7 — 

£ 0 +4£i ] 

1 c 2j 

_ 

q 

.v c l, 

. 

c 0 

V C 0J 

JJ 


(33a) 


(33b) 


n yn (33c) 
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y 4 


(Air 


+- 


i 

c 3 

1 

c 0 


3-—) 

V c 3) 

( 


£ 3+64 


+ - 


C2 


21-10 


V 


q_ 
c o ) 


£q -1-lOq 


' 9-4 s' 

£ 2 +4£ 3 

1 

+ — 

Tl5-7^) 

£l + 7£ 2 

l c 2) 
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a q j 




(33d) 



( Al ) 2 
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1 

+ 
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+ — 
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3-^ 

£4 +£ 5 

9-4 — 

£ 3 +4£ 4 

15-7 — 

£2 + 7£ 3 
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.V C 4 , 


c 3 

A c 3 ) 


c 2 

A c 2 / 



Y ■) 

21 — 10 — 

+ 10^2 

1 

+ — 

Y ■) 

27-13— 

£q+13£i ] 

A c i ) 


c 0 

A c oJ 

JJ 


(33e) 


B (4M> 

y n = 



3(2y — 1) — (3 j — 2) - 


c n— j+l 
c n—j 


^ n—j "I" (3j 1-1 [ n ^ n ~^ 


(33f) 
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If the right end is also simply supported, the right-end slope, tan 6 n , at the strain-sensing station, x n , is 
nonzero ( tan 0„ ^ 0 ), and slope equation (19) can be used to determine the unknown right-end slope, 
tan , in terms of the left-end slope, tan 0 q , calculated from equation (32) as 


„ 4/ V 1 

tan0 n = — Zj — 


2 j=l c i~' 

n 


2 - 


c j 

C M 


£j-l+£; 


+ tan 0 q 


_ Al A 1 


2 - 


c j 


C J- 1 ; 


£j-l+£j 


nAl 


y n 


(34) 


Note that the right-end slope, tan0 rt , does not appear in deflection equation (33); therefore, the knowledge 
of tan 0 H is not needed in the deflection calculations. 


Left End Fixed 

Deflection equation (23) for the nonuniform cantilever beam is also applicable to the nonuniform simple 
beam with the left end fixed (the right end can be either fixed or simply supported). The deflection, y n , at 
the right end, calculated from cantilever beam deflection equation (23 f) with exact strain inputs, 
theoretically should be zero. If the calculated right-end deflection, y n , should show a small nonzero error, e 
(y„ =e«l), however, then deflection equation (23) for the cantilever beam ( yo = tan0q = 0) must be 
corrected for the present case as 


B W 

yi =■ 


1 i 

i- 


3(2 j — 1) — (3j — 2) 


c i~j + 1 

C i~J 


£j—j+(3j 2 )£j_j + \ 


yj [ equation (23); yo=tan0o=O] 


i 

— e 

Correction 

term 


(35) 


Similar to figure 4, equation (35) is used to proportionally shift the calculated deflection curve slightly 

D 

downward or upward depending on the sign of e to make the right-end deflection zero (y„ = 0). Keep in 
mind that the degree of shifting for the present case is minuscule compared with the previous case (left end 
simply supported, eq. (33)). 

TWO-LINE STRAIN-SENSING SYSTEM 

The two-line strain-sensing system proposed by Ko (ref. 2) is a powerful method for simultaneously 
monitoring the bending and torsional deformation of nonuniform cantilever beams under combined bending 
and torsion. The two-line strain-sensing system eliminates the need for installing the shear strain sensors to 
measure the surface distortions through which the cross-sectional rotations of the beam can be determined. 

Sensor Installation Examples 

Installation of the two-line strain-sensing system for typical beam-like structures is discussed in the 
succeeding sections. These typical beam-like structures include a wing box and cantilever tubular beam. 
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Wing Box 

Figure 5 shows the first example of applying the two-line strain-sensing system to a wing box structure 
under combined bending and torsion. The two strain-sensing lines (front and rear) can be installed along the 
lower edges (front and rear) of the wing box with chord-wise separation distance, dj , at the strain-sensing 
cross section, x = Xj . The two strain-sensing lines can be parallel ( dj = constant) or tapered with the 
separation distance, dj , decreasing toward the wingtip according to the horizontal taper rate of the wing. 

Cantilever Tubular Beam 

Figure 6 shows another example of the application of the two-line strain-sensing system to a tapered 
cantilever tubular beam under bending and torsion. The front and rear strain-sensing lines can be installed at 
a symmetrical angular location of ±a measured from the bottom of the tube. To obtain larger bending strain 
data, a smaller value of a is desired. But to obtain more accurate cross-sectional twist data, a larger a (that 
is, a larger separation distance) is desired. Therefore, a = ±30 ° may be a good compromise, because the 
separation distance, dj , between the two strain-sensing lines will be equal to the tube radius. The magnitude 
of the strains measured at a = ±30 ° location (numerators of equation (23)) will be cos a times the strains 
(numerators of equation (23)) measured at the bottom of the beam. The vertical distance, Cj, from the 
neutral axis to the strain-sensing stations (denominators of equation (23)), however, will be cos a times the 
tube radius, thereby canceling out the effect of cos a . Therefore, the averaged deflection curve generated by 
using the two-line strain-sensing system at a = ±30 ° will be equivalent to the deflection curve generated by 
using a single strain-sensing line along the bottom surface of the tube. 

Deflections 

The deflections, { >7 , y \ } , of the front and rear strain-sensing stations lying at the sensing cross section, 
x = Xj (i = 0, 1, 2, 3, ..., n) (figs. 5, 6), can be calculated from deflection equation (23), which is rewritten 
(terms vq + (i)Al tan (9 q removed) in the following forms for the front and rear strain-sensing stations, 
respectively: 

For the front strain-sensing line, 


y, = 


(AD 


2 i , 

-Y — 

6 


3(2 j -1)- (3y — 2)- 


c i~j+ 1 
C i~j 


£ i-j + (3j -2)£j-j + if ; (i = 1, 2, 3, ..., ri) (36) 


For the rear strain-sensing line, 


, AD 2 ± 1 

»■ = — 7 — L 

6 j=i c i~j 


3(2; -1)- (3;- 2)- 


c i~j+ 1 
c i~j 


£ i-j + (3 j ~ 2)£'j_ j+l ; (* = 1,2,3 n) (37) 


By averaging the front and rear deflections, { yj,y\ }, calculated from equations (36) and (37), respectively, 
we can eliminate the torsional effect to yield the true bending deflections of the wing box. 
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Cross-Sectional Twists 


If dj denotes the chord-wise separation distance between the front and rear strain-sensing cross section 
at x =X[ (figs. 5, 6), then the cross-sectional twist angle, 0; , at the strain-sensing cross section, x =Xj , can 
be calculated from the cross-sectional twist angle equation 


0/ = sin 




yj-yt 

c d-i j 


(i = 0, 1,2,3, ...,«) 


(38) 


in which 0 q = 0 at the wing root. The deflections, { y\ ,y'j}, can be calculated from deflection equations 
(36) and (37), respectively, using the measured bending strain data as inputs. 

For an aircraft wing (fig. 5), when the value of dj is much larger than the wing depth, equation (38) 
could give sufficiently accurate 0/ without the need to install distortion sensors (see the section entitled, 
“Tapered Wing Boxes”), because the effect of surface shear strains on the bending strains is negligible. For 
the tubular structure (fig. 6) for which dj is in the order of the radius, however, equation (38) may lose 
accuracy, because the surface shear strains will influence the bending strains (see the section entitled, 
“Tapered Tubular Beams”). An alternative method to obtain the cross-sectional twist angles, 0,- , for a 
tubular structure is to install a surface shear strain-sensing line for measurement of surface distortion angles, 
from which 0/ can be calculated. This alternative method is discussed in reference 2. 

APPLICATION TO FREE-FREE BEAMS 

The fuselage of an aircraft during flight is a good example of a free-free unsupported beam (fig. 7). The 
fuselage is a tubular structure with varying cross sections and can be considered as a nonuniform beam so 
that deflection equations (16) and (23), developed for nonuniform and slightly nonuniform cantilever beams, 
respectively, can be applied. For monitoring the deflections of the airborne fuselage under either vertical or 
horizontal bending, the strain-sensing line can be installed axially on the bottom (belly) surface or on the 
side of the fuselage (fig. 7). 

Before deflection equation (16) or (23) is applied for the deformed shape predictions of the fuselage, 
some basics of equations (16) and (23) must be understood Remember that in the formulation of deflection 
equations (16) and (23), the reference strain-sensing station, xq, is located at the built-in end. With the 
built-in condition, yo = tan0o=O, imposed, equations (16) and (23) are also applicable to the moving 
cantilever beams, because the deflection, yj , calculated from equation (23) is relative to the tangent line 
( tan 0() = 0) stemming from the reference strain-sensing station, xq, at the built-in end (fig. 7). Equations 
(16) and (23) can be applied to the free-free unsupported moving fuselage case using the following 
approach. 

The fuselage can be considered two nonuniform cantilever beams joined together at the center-of- 
gravity (CG) cross section (fig. 7). During flight, the CG cross section has the least movement compared 
with other fuselage cross sections. For practical purpose, it is best to choose the reference strain-sensing 
station, xq, at the CG cross section and divide the strain-sensing line into two segments (fig. 7). One 
segment is for the fore section and the other for the aft section of the fuselage. Thus, by imposing 
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vq = tan Oq = 0 at the reference strain-sensing station, xq (considered as a built-in end), deflection 
equations (16) and (23) for the cantilever beam can be used to calculate the deflections of the fore and aft 
sections of the fuselage. The calculated deflection, y ; - , will then be the relative deflection with respect to the 
tangent line passing though the reference strain-sensing station, xq , which moves (translates and rotates) 
with the fuselage (fig. 7). 

APPLICATION TO STRUCTURAL PANELS 

Deflection equations (33) and (35) developed for slightly nonuniform simple beams can also be used for 
deformed shape predictions of structural panels (for example, flat, cylindrical, shallow spherical, and wing 
panels). For the shape sensing of structural panels, multiple biaxial strain-sensing lines are required (fig. 8). 
The panel strip region lying beneath each strain-sensing line will behave as an embedded beam (cantilever 
or simple beam, depending on the panel edge support conditions); therefore, deflection equations (33) and 
(35) are applicable. Each surface strain sensor will then sense the bending strain of the fictitious beam 
embedded beneath the strain-sensing line and not sense the bending strain of an isolated beam. As such, 
each strain output is the true panel bending strain with a biaxial effect. Keep in mind that the bending strain 
of the beam that is embedded in the panel will increase if the embedded beam is isolated. Therefore, the 
one-dimensional beam strain does not need to be corrected into a two-dimensional panel strain. This feature 
is a powerful characteristic of the Ko displacement theory, which is applicable whether the structure is one- 
dimensional and beam-like or two-dimensional and panel-like. The deflection equations previously 
developed for uniform beams (cantilever and two-point supported; ref. 2) were used for the deformed shape 
predictions of uniform-thickness square panels (four edges simply supported or four edges clamped), and 
were quite accurate in the panel shape predictions. 

VALIDATION OF KO DISPLACEMENT THEORY 

To examine the accuracy of the further-developed Ko displacement theory for nonuniform beams, 
several beam structure cases were chosen: (1) tapered tubular beams (cantilever and two-point supported), 
and (2) wing boxes (depth-tapered unswept, depth-tapered swept, width-tapered unswept, and double- 
tapered unswept). The Structural Performance And Resizing (SPAR) finite-element computer program 
(ref. 7) was then used to generate both the nodal strains and nodal displacements for these example 
structures. All the structures considered were made of aluminum material. 

Input Strains 

To test the accuracy of the displacement theory for nonuniform beams, the fiber-optic measured strain 
data were needed for input to the deflection equations for the calculations of beam deflections. Because the 
measured bending strain data were not yet available from the experiments, the SPAR program was used to 
generate the bending strains. The SPAR bending strains were generated as follows. 

If (7 j is the axial stress at the strain-sensing station, x ( - , obtained from the SPAR nodal stress outputs, 
then the corresponding bending strain, £*• , for input to the deflection equations can be calculated from 
Hooke’s law. 


e; = 


E 


(39) 


20 



Using the finite-element generated strains (eq. (39)), we can calculate the structural deflections and then 
compare them with the corresponding nodal displacements obtained from the SPAR displacement outputs 
for validating the accuracy of the Ko displacement theory. 

Alternatively, using the SPAR displacement outputs, we can also obtain the bending strain, £ L , by 
dividing the element axial length change by the undeformed length of the element at the strain-sensing 
station, X( (ref. 2). This element displacement method gives bending strains that are practically identical to 
those obtained from the nodal stress method (eq. (39)) in the region of small deflections, but it starts to lose 
accuracy in the highly bent region near the free end. The SPAR displacement outputs are referred to the 
undeformed beam axis, and at increasing deflections near the highly bent beam tip, the true tensile strains 
(positive value) along the deformed beam axis, when projected on the undeformed beam axis, may appear as 
compressive strains (negative value). Thus, in the calculations of strains using the element displacement 
method, y-displacement components must be considered. Therefore, the simpler nodal stress method 
(eq. (39)) was used in generating the input strains. 

Tapered Tubular Beams 

Each of the tapered tubular beams (cantilever, two-point supported) used in the SPAR analysis has a 
radius that linearly decreases from one end toward the other end. The tapered tubular beam tapers down not 
only in the vertical and horizontal planes, but also in any radial plane because of circular cross sections. 
Table 1 lists the dimensions of each aluminum tapered tubular beam considered. 

Table 1. Dimensions of aluminum tapered tubular beams. 


/, in. 

co , in. 

c n , in. 

c n /co 

100.5 

4 

4 

1 .00 (uniform) 

100.5 

4 

3 

0.75 

100.5 

4 

2 

0.50 

100.5 

4 

1 

0.25 


Cantilever Beams 

Figure 9 shows the finite-element model for a tapered cantilever tubular beam. The beam is subjected to 
a tip vertical load of P = 100 lb (fig. 9a), and a combined tip vertical load and tip torque of P = 100 lb and 
T= 400 xc„ in-lb, respectively (fig. 9b). 

Bending Only 

Figure 10 shows the strain curves for different tapered cantilever tubular beams under bending only, 
calculated from equation (39), using the SPAR element stress outputs. The strain curve is a straight line for 
the classical uniform beam case (c n / cq = 4/4) but becomes bow shaped as the depth ratio, c n / cq > 
decreases. When the depth ratio reaches a minimum of c n I cq =1/4, the strain curve is deeply bent by 
nearly 100°. 

Figure 11 shows the beam deflection curves for different tapered cantilever tubular beams, calculated 
from deflection equation (16) (eq. (23) for the uniform beam), using the strain data at the strain-sensing 
stations shown in figure 10 (solid curves with solid circular symbols). For comparison, the deflections 
calculated from the SPAR program are also plotted (dashed curves with open circular symbols). The 
excellent correlation between the two sets of deflection curves shows the high degree of accuracy of the 
displacement theory (eqs. (16) and (23)) formulated for the nonuniform cantilever beams. 
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Keep in mind that in the formulation of the Ko displacement theory (ref. 2), the distribution of bending 
strains along the discretized cantilever beam was assumed to be piecewise linear as described by 
equation (9). Therefore, the strain curves associated with this theory are the piecewise linear approximations 
(connecting the strain data points at the strain-sensing stations with straight lines) of the smooth strain 
curves shown in figure 10. Note from figure 10 that for slightly bent strain curves, the piecewise linear 
approximation could be quite accurate, but in a highly bent region of the strain curve for the lowest depth 
ratio, c n / cq =1/4 , the piecewise linear approximation may lose some accuracy locally. In spite of this 
shortcoming, the Ko displacement theory (eqs. (16) and (23)) still predicts highly accurate deflections for 
the whole range of depth ratios, 4 /4 > c n Icq > 1 /4 (fig. 11). The same arguments also hold for the 
remaining examples presented in the succeeding sections. 

Bending and Torsion 

Figure 12 shows the SPAR-generated strain curves for the front (fig. 12a) and rear (fig. 12b) strain- 
sensing lines for different tapered cantilever tubular beams under combined bending and torsion. The front 
(left) and rear (right) strain-sensing lines are 60° apart as shown in figure 6. The front and rear strain curves 
are very similar because of the proximity of the two strain-sensing lines. Note from figure 12 that the strain 
curves for the tapered tubes (c n / cq < 4/4) are bow shaped, and the degree of bend increases as the depth 
ratio, c n / cq , decreases. When the depth ratio reaches a minimum of c n I cq =1 / 4 , the front and rear strain 
curves are deeply bent by nearly 100°, similar to the strain curves in the bending only cases (fig. 10). For the 
classical uniform cantilever beam case (c n / cq = 4 / 4 , fig. 12), the front and rear bending strain curves are 
linear except in the beam tip region. The abnormal deviation of bending strains (especially for the front 
strain curve) near the beam tip region could be attributed to the increasing local element distortions caused 
by the concentrated tangential loads applied at selected points at the beam tip to create torsion. Keep in mind 
that the bending strains have been calculated from the finite-element stresses. 

If torsion is removed, the front and rear strain curves of the classical uniform cantilever beam case 
(c n / cq = 4/4) will be perfectly straight lines as shown with dashed lines in figure 12 (see also fig. 10 for 
the c n Icq =4/4 case). With torsion added, each bending strain curve for the c n t cq =4/4 case (fig. 12) 
gradually deviates from the corresponding strain curve for bending only (dashed lines, fig. 12), and the 
deviation increases slightly toward the beam tip. 

Figure 13 shows the deflection curves along the front and rear strain-sensing lines for different tapered 
cantilever tubular beams, calculated from deflection equation (16) (eq. (23) for the uniform beam), using the 
strain data shown in figure 12 (solid curves with solid circular symbols). For comparison, the deflections 
calculated from the SPAR program are also plotted (dashed curves with open circular symbols. The 
excellent correlation between the two sets of deflection curves demonstrates the high accuracy of the 
displacement theory (eqs. (16) and (23)) formulated for nonuniform cantilever beams. 

Figure 14 shows the cross-sectional twist angles, 0; , for different tapered cantilever tubular beams, 
calculated from equation (38) using the deflection data, { v* , y \ } , along the front and rear strain-sensing 

lines shown in figure 13. Graphically, the agreement between the twist angles calculated from the 
displacement theory and those calculated from the SPAR program does not appear as good as the agreement 
between the deflection curves (fig. (13)). It must be understood that the values of 0; are quite small (less 

than 1°) because of the proximity (1-4 in. apart) of the front and rear strain-sensing lines, and because of a 
low twisting moment, T = 400 xc„ in-lb. Therefore, slight differences (almost graphically inconspicuous) 
between the Ko displacement theory and SPAR program deflection data of figure (13) have a strong effect 
on the values of 0/ calculated from equation (38). Table 2 lists the beam tip deflections, { y n , y' n }, of the 
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front and rear strain-sensing lines, and the cross-sectional rotation angles, 0; , calculated from the SPAR 
program and Ko displacement theory. 

Table 2. Comparison of displacement data calculated from SPAR program with those calculated 
from Ko displacement theory for different tapered cantilever tubular beams; / = 100.5 in., P = 1 00 lb, 
T ~ 400 X c n in-lb. 


CO .in. 

c n , in. 

Jk- 

in. 

y'n . 

in. 

0/7 5 

deg 

SPAR 

Ko 

SPAR 

Ko 

SPAR 

Ko 

4.0* 

4.0* 

0.72027 

0.72677 

0.70249 

0.71001 

0.25468 

0.24007 

4.0 

3.0 

0.88400 

0.88706 

0.86854 

0.87291 

0.29526 

0.27025 

4.0 

2.0 

1.16457 

1.16199 

1.15133 

1.15148 

0.37930 

0.30109 

4.0 

1.0 

1.78962 

1.77362 

1.77859 

1.76628 

0.63197 

0.42056 


* Uniform beam. 


Note from table 2 that the differences between the <j) n values calculated from the SPAR program and 
those calculated from the Ko theory are minuscule, in the range of {0.0146, 0.0250, 0.0782, 0.2114} deg, 
for the depth ratios, c n /cq = {4 / 4, 3/4, 2/4, 1 / 4} , respectively. As seen in the subsequent wing box cases, 

in which the values of 0/ are larger because of greater distances between the front and rear strain-sensing 
lines, the prediction of 0/ turns out to be quite good. 

Two-Point Supported Beams 

Two cases of simple beams (two-point supported beams) were considered. Figure 15 shows the SPAR 
models of the two cases of tapered tubular simple beams, one with both ends simply supported (fig. 15a), 
and one with both ends fixed (fig. 15b). Each simple beam is subjected to a downward load of P = -200 lb 
at the center of the beam. 

Simply Supported Ends 

Figure 16 shows the strain curves, calculated from equation (39) using the SPAR element stress outputs, 
for different tapered tubular simply supported beams subjected to a downward force of P = -200 lb at the 
center of the beam. The strain curve is roof shaped for the classical uniform beam (c n I cq = 4 / 4 ). For the 
depth ratio, c n / cq < 4 / 4 , the strain curves are slightly concave in the left-hand region but convex (bow 
shaped) in the right-hand region with the degree of bow bend increasing with the decreasing depth ratio, 
c n Icq. At the lowest depth ratio of c n / cq = 1 / 4 , the strain curve in the right-hand region (decreasing 
radius) is deeply bent by nearly 90°. 

Figure 17 shows the beam deflection curves for different tapered tubular simply supported beams, 
calculated from simply supported beam deflection equation (33) using the strain data shown in figure 16 
(solid curves with solid circular symbols). These deflection curves are compared with the deflection curves 
calculated from the SPAR program (dotted curves with open circular symbols). The SPAR program gives 
slightly larger deflections, because the transverse shear effect is accounted for in the SPAR program. Note 
from figure 1 7 that the transverse shear effect (difference between the two sets of deflection curves) is quite 
small for the uniform beam, c n Icq = 4 14 , and gradually grows larger as the depth ratio, c n Icq , decreases. 
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At the lowest depth ratio, c n / eg = 1 / 4 , the transverse shear effect at the center of the beam amounts to 
3.63 percent of the beam center deflection. Remember, for the tapered cantilever tubular beam cases 
(fig. 11), the transverse shear is practically inconspicuous for the whole range of depth ratios, 
4 / 4> c n / cq >1 / 4 . 

To show that the difference between the deflections calculated from the SPAR program and those 
calculated from the displacement theory is caused by the transverse shear effect, the shear modulus, G, was 
set very large ( G — > « ) in the SPAR program to eliminate the transverse shear effect. As expected, the 
resulting SPAR deflection curves (dashed curves with open square symbols, fig. 17) practically fall on top 
of the corresponding Ko deflection curves (solid curves with solid circular symbols, fig. 17), thus proving 
the argument. The excellent deflection predictions based on deflection equation (33) for the simply 
supported beam validate the accuracy of the Ko displacement theory developed for slightly nonuniform 
beams. 

Fixed Ends 

Figure 18 shows the strain curves generated from the SPAR program (equation (39)), for different 
tapered tubular simple beams with two ends fixed, subjected to a central downward load of 
P = -200 lb. The strain curve is again roof shaped for the classical uniform beam (c n Icq = 4 /4 ). For the 
depth ratio, c n / cq < 4 / 4 the strain curves are slightly concave in the left-hand region but convex in the 
right-hand region, with the degree of convexity increasing as the depth ratio, c n / cq , decreases. At the 
lowest depth ratio of cq / c n = 1 / 4 , the strain curve in the right-hand region is deeply bent. 

Figure 19 shows the beam deflection curves for the fixed-end simple beams, calculated from cantilever 
beam deflection equation (16) (eq. (23) for the uniform beam), using the strain data shown in figure 18 
(solid curves with solid circular symbols). These deflection curves are compared with the deflection curves 
calculated from the SPAR program (dotted curves with open circular symbols). Note that equations (16) and 
(23) give y n = 0 at the right end, and therefore it is not necessary to use the error correction equation (35) 
for the y M ^0 case. Note from figure 19 that the correlation between the two sets of deflection curves is 
reasonably good. Similar to the simply supported case, the SPAR program gives slightly larger deflections, 
because the transverse shear effect is considered in the SPAR program. Note that the transverse shear effect 
(difference between each set of SPAR and Ko deflection curves) increases slightly as the depth ratio, 
c n Icq , decreases. At the lowest depth ratio, c n / cq = 1/4, the transverse shear effect at the center of the 
beam amounts to 10.29 percent of the beam center deflection. Again, by setting the shear modulus, G, very 
large ( G — > °° ) in the SPAR program, the transverse shear effect was removed, and the resulting SPAR 
deflection curves (dashed curves with open square symbols) fall nearly on top of the corresponding Ko 
deflection curves (solid curves with solid circular symbols) for the whole range of depth ratios, 
4/4 >c n /cQ > 1 / 4 . Thus, if the transverse shear effect is ignored, then cantilever beam deflection 
equations (16) and (23) are quite accurate when applied to the fixed-end tapered simple beam, even down to 
the lowest depth ratio of cq / c n = 1 / 4 . 


Tapered Wing Boxes 

For numerical validation of Ko displacement theory, three different types of wing boxes were 
considered: (1) depth-tapered wing boxes (unswept and swept), (2) width-tapered unswept wing boxes, and 
(3) double-tapered unswept wing boxes. The two strain-sensing lines for each wing box case are installed at 
the front and rear lower edges of the wing box. 
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Depth-Tapered Wing Boxes 


The depth-tapered wing boxes considered are unswept (fig. 20a) and swept (swept angle, /), deg; 
fig. 20b). All depth-tapered wing boxes have a constant width, with the depth tapering off toward the free 
end (vertically tapered). Table 3 lists the dimensions of the depth-tapered unswept and swept wing boxes. 
Each wing box is subjected to vertical load of P = 100 lb and a torque of T = 50 xw in-lb (clockwise) at the 
tip of wing box (fig. 20). For the swept wing boxes, a counterclockwise torsion of T = -50 xw in-lb was 
also considered. 


Table 3. Dimensions of depth-tapered unswept and swept wing boxes. 


l, in. 

c 0 , in. 

c n , in. 

w , in. 

A de g 

(unswept) 

(swept) 

100.0 

2.0 

2.0 (uniform) 

17.0 

0.0 

10.0 

100.0 

2.0 

1.5 

17.0 

0.0 

10.0 

100.0 

2.0 

1.0 

17.0 

0.0 

10.0 

100.0 

2.0 

0.5 

17.0 

0.0 

10.0 


Unswept Wing Boxes 

Figure 21 shows the plots of the strains along the front and rear strain-sensing lines, calculated from the 
SPAR element stress outputs, for different depth-tapered unswept wing boxes subjected to a vertical load of 
P = 100 lb and a torque of T= 50 xw in-lb at the wing box tip (fig. 20). The front-edge and rear-edge strain 
curves for the nontapered uniform wing box ( c n / cq = 2 / 2 ; fig. 2 1 ) are practically straight lines like those 
in the classical cantilever beam case. For the tapered wing boxes, the front-edge strain curves 
(fig. 2 1 a) become bow shaped with the bend increasing as the taper depth ratio, c n / cq , decreases (very 
similar to the tapered cantilever beam cases, fig. 10). When the taper depth ratio reaches a minimum of 
c n / c () = 0.5/2, the strain curve becomes deeply bent by more than 90°. The rear-edge strain curves for the 
tapered wing boxes are almost straight except in the wing root and wingtip regions. As the taper depth ratio, 
c n / cq , decreases, the strain curves become more tilted. 

Figure 22 shows the front-edge and rear-edge deflections, { >’/ , yj } , for different depth-tapered unswept 

wing boxes, calculated from deflection equations (36) and (37), respectively, using the strain data shown in 
figure 21 as inputs (solid lines with solid circular symbols). For comparison, the deflections obtained from 
the SPAR nodal displacement outputs are also plotted (dashed lines with open circular symbols). The 
excellent correlation between the two sets of deflection curves reinforces the high degree of accuracy of the 
Ko displacement theory (eqs. (16) and (23)) formulated for nonuniform cantilever beams. 

Note that under the combined bending and torsional loading, the front-edge and rear-edge deflection 
curves shown in figure 22 contain torsional components. To obtain an overall bending deflection curve for 
the wing box, the torsional effect must be removed by averaging the front-edge and rear-edge deflection 
curves of figure 22. 

Figure 23 shows the cross-sectional twist angles, (f>i , for different depth-tapered unswept wing boxes, 
calculated from equation (38) using the deflection data, { y,- , y\ } , along the front and rear strain-sensing 
lines. Again the agreement between the displacement theory and SPAR program is quite good. Interestingly, 
at the taper depth ratios of c n Icq = 0.5/2 and c n Icq = 1/2, the cross-sectional twist angle, (under 
combined loading), in the inboard region has the opposite sign of the wingtip torque, T. This sign change 
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of 0/ could result from the lack of cross-sectional shear stiffness, because the SPAR wing box models are 
hollow with no interior shear reinforcements. It was found that if rib and spar reinforcement elements were 
added to the hollow SPAR models, the cross-sectional twist angle, 0/ , at any point along the wingspan, 
would follow the sign of the wingtip torque, T. 

Swept Wing Boxes 

Figure 24a shows the plots of the strains along the front strain-sensing line, calculated from the SPAR 
element stress outputs, for different depth-tapered swept wing boxes subjected to a vertical load of 
P = 100 lb and a torque of T = 50 X w n in-lb (clockwise, looking toward the wing root) at the wing box tip 
(fig. 20b). Overall the shapes of the front-edge strain curves are similar to those for the unswept wing boxes 
cases (fig. 21) except in the wing root regions where the strain curves have saddle-shaped areas. The 
front-edge strain curve for the uniform wing box ( c n / cq = 2 / 2 ; fig. 24a) is a straight line except in the 
wing root region. For the depth-tapered swept wing boxes, the front-edge strain curves become bow shaped 
(except in the wing root region), with the bend increasing as the taper depth ratio, c n / cq , decreases (very 
similar to the depth-tapered unswept wing box cases, fig. 21). When the taper depth ratio reaches a 
minimum of c n / cq = 0.5/2, the strain curve becomes deeply bent by more than 90°. 

Figure 24b shows the plots of the strains along the rear-edge strain-sensing line, calculated from the 
SPAR element stress outputs, for different depth-tapered swept wing boxes. Overall the rear-edge strain 
curves are almost linear like those for the unswept wing boxes cases (fig. 21) except in the wing root regions 
where the strains increase sharply. The rear-edge strain curve for the uniform wing box ( c n / cq = 2 / 2 ) is 
almost a straight line except in the wing root and wingtip regions. For the tapered swept wing boxes, the 
strain curves are nearly straight except in the wing root and wingtip regions. As the taper depth ratio, 
c n / cq , decreases, the strain curves become more tilted. 

Figure 25 shows the front-edge and rear-edge deflections, { }, for different depth-tapered swept 

wing boxes under bending and clockwise torsion, calculated from deflection equations (36) and (37), using 
the strain data shown in figure 24 as inputs (solid lines with solid circular symbols). For comparison, the 
deflections obtained from the SPAR nodal displacement outputs are also plotted (dashed lines with open 
circular symbols). Note that the correlation between the Ko and SPAR deflection curves is excellent, 
providing confidence in the Ko displacement theory (eqs. (16) and (23)) formulated for nonuniform 
cantilever beams. 

Figure 26 shows the cross-sectional twist angles, 0/ (with respect to airstream), for different 
depth-tapered swept wing boxes (clockwise torsion cases), calculated from equation (38) using the 
deflection data, { yf , y \ } , along the front-edge and rear-edge strain-sensing lines. Similar to the unswept 

wing cases (fig. 23), at the taper depth ratios, c n / cq = 0.5/2 and c n / cq = 1/2, the cross-sectional twist 
angle, 0/ , in the inboard region has the opposite sign of the wingtip torque, T, because the SPAR wing box 
models are hollow, having no interior reinforcements (no cross-sectional shear stiffness). The peculiar 
behavior of 0/ (under combined loading) disappeared if the interior shear reinforcement elements were 
added in the hollow SPAR models. Overall the agreement between the Ko displacement theory and the 
SPAR program is excellent for the swept wing boxes. 

Figures 27, 28, and 29 show sets of strain curves, deflection curves, and cross-sectional twist angle 
curves, respectively, for the swept wing boxes under combined bending and counterclockwise torsion. 
Notice that by reversing the direction of torsion, the shapes of the front and rear strain and deflection curves 
do not become the interchanged versions of the corresponding curves for the clockwise torsion cases 
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(figs. 24, 25, 26) but have slight variations resulting from the swept effect near the wing root and wingtip. 
The most conspicuous difference between the two torsion cases is that the front-edge strain curve of 
c„ = 0.5 in. for the counterclockwise torsion case (fig. 27b) does not have a saddle-shaped area near the 
wing root like the rear-edge strain curve of c n = 0.5 in. for the clockwise torsion case (fig. 24a). The 
excellent agreement shown in figures 28 and 29 reinforces the high degree of accuracy of the Ko 
displacement theory (eqs. (16) and (23)) formulated for nonuniform cantilever beams. 

Width-Tapered Wing Boxes 

Figure 30 shows the unswept width-tapered wing box with a constant half depth, cq = c n = 0.5 in., and a 
decreasing width toward the free end (horizontally tapered). The strain-sensing lines are located at the lower 
front and rear edges. Each beam is subjected to a vertical load of P = 100 lb and a torque of T = 50 x w n 
in-lb (clockwise) at the wing box tip (fig. 30). Table 4 lists the dimensions of each width-tapered unswept 
wing box. 


Table 4. Dimensions of width-tapered unswept wing boxes. 


l, in. 

c 0 , in. 

c n , in. 

wq , in. 

w n . in - 

/3,deg 

60.0 

0.5 

0.5 

30.0 

30.0 ( uniform ) 

0.0 

60.0 

0.5 

0.5 

30.0 

20.0 

0.0 

60.0 

0.5 

0.5 

30.0 

10.0 

0.0 

61.0 

0.5 

0.5 

30.0 

5.0 

0.0 


Figure 31 shows the strain curves generated from the SPAR element stress outputs for different width- 
tapered wing boxes subjected to a vertical load of P = 100 lb and a torque of T = 50 x w n in-lb (clockwise) 
at the wing box tip (fig 30). 

Figure 32 shows the deflection curves for different width-tapered unswept wing boxes, calculated from 
deflection equation (23), using the associated input strain data of figure 31 (solid lines with solid circular 
symbols). For comparison, the deflection curves obtained from the SPAR nodal displacement outputs are 
also plotted (dashed lines with open circular symbols). The excellent correlation between the two sets of 
deflection curves validates the high degree of accuracy of the Ko displacement theory (eq. (23)) formulated 
for nonuniform cantilever beams. Figure 33 compares the cross-sectional twist angles, </>/ , calculated from 
equation (38) with those generated from the SPAR program, for different width-tapered wing boxes. Again 
the agreement is excellent. 

Double-Tapered Wing Boxes 

The double-tapered wing box (fig. 34) is unswept and has both a width and depth (horizontally and 
vertically) that taper down toward the free end. The front and rear strain-sensing lines are located at the front 
and rear lower edges, respectively, of the wing box. Table 5 lists the dimensions of the double-tapered wing 
box. The wing box is subjected to a vertical load of P = 100 lb and a torque of T = 50 xw n in-lb (clockwise) 
at the wing box tip. 
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Table 5 . Dimensions of double-tapered unswept wing box. 


l, in. 

c 0 , in. 

c n , in. 

wq , in. 

w n , in. 

P • deg 

100.0 

2.0 

1.0 

34.0 

17.0 

0.0 


Figure 35 shows the strain curves generated from the SPAR element stress outputs for the front and rear 
oblique strain-sensing lines at the lower edges of the double-tapered wing box. The strain curves are deeply 
bent because of the double tapering effect. Remember that for the classical uniform cantilever beam, the 
strain curve is a straight line. 

Figure 36 shows the deflection curves along the two oblique strain-sensing lines, calculated from 
deflection equation (16), using the strain data generated from the SPAR program shown in figure 35 (solid 
lines with solid circular symbols). For comparison, the deflections calculated from the SPAR program are 
also plotted (dashed lines with open circular symbols). The excellent correlation between the two sets of 
deflection curves validates the high degree of accuracy of the Ko displacement theory (eq. (16)) formulated 
for nonuniform cantilever beams. 

Figure 37 compares the cross-sectional twist angles, fa , calculated from equation (38) with those 
calculated from the SPAR program, for the double-tapered wing box. The Ko and SPAR curves for fa 
practically fall on top of each other, showing the surprisingly high degree of accuracy of the Ko 
displacement theory. 

EXACT AND EXPANDED DEFLECTION EQUATIONS 

All the deflection curves for vertically tapered beams predicted from the Ko displacement theory, as 
shown in the preceding sections, were calculated using the exact deflection equation (16) formulated for 
nonlinear beams (for example, tapered beams). Because deflection equation (16) breaks down for the 
uniform beams, the deflection curves for the uniform beam cases had to be calculated from the expanded 
deflection equation (23), which was obtained from the exact deflection equation (16) by expanding the 
logarithmic term, logqTq / q _ | ) , in the neighborhood of (q / q_\)~ 1 . The preceding examples demonstrate 
that equation (16) provides quite accurate shape predictions for tapered beams with all the taper depth ratios 
down to c n / cq =1 / 4 . 

To examine the deflection prediction accuracy, expanded deflection equation (23) was used to calculate 
the deflections of tapered cantilever tubular beams with different taper depth ratios, as shown in figure 11. 
Table 6 compares the beam tip deflections calculated from the SPAR program with those calculated from 
the exact deflection equation (16) and those calculated from the expanded deflection equation (23). 

Note from table 6 that the two sets of deflection data calculated from deflection equations (16) and (23) 
are very close, with a peak difference of only 0.5293 percent for the lowest depth ratio, c n / cq = 1/4 . Thus, 
if the deflection data of equation (23) are plotted out in figure 11, the deflection curve will fall practically on 
top of the deflection curve calculated from equation (16) (fig. 11). This observation indicates that the 
expanded deflection equation (23) is quite accurate for the whole range of depth ratios, 
4 / 4 > c n / cq>\! 4 . 
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Table 6. Comparison of beam tip deflections, y n , of tapered cantilever tubular beams calculated from SPAR 
program with those calculated from displacement equations (16) and (23); l =100.5 in., P =100 lb, n = 8. 


co , in. 

c n , in. 

SPAR 

y n , in. 
Eq. (16) 

(exact) 

Eq. (23) 

(expanded) 

Percent difference 
between 

eq. (16) and eq. (23) 

4.0 

4.0 (uniform) 

0.711620 

Not applicable 

0.702304 


4.0 

3.0 

0.876340 

0.865681 

0.866592 

0.1052 

4.0 

2.0 

1.158010 

1.145292 

1.145487 

0.0170 

4.0 

1.0 

1.784170 

1.762277 

1.752950 

-0.5293 


CONCLUDING REMARKS 

The Ko displacement theory, originally formulated (in recursion formats) for shape predictions of 
nonuniform cantilever beams, was further developed into more compact mathematical forms. The accuracy 
of the further-developed Ko displacement theory was validated by finite-element analyses of various beam 
structures. The principal findings are as follows: 

1. In the further-developed Ko displacement theory for the nonuniform cantilever beam, the 
deflections at the strain-sensing stations were expressed explicitly in terms of geometrical 
parameters of the beam and bending strains at all the inboard strain-sensing stations. 

2. The exact deflection equations for nonuniform cantilever beams were expanded for slightly 
nonuniform cantilever beams so that the expanded deflection equations could degenerate into the 
limit case for uniform cantilever beams without encountering mathematical breakdown problems. 

3. For tapered cantilever tubular beams with different taper rates, the deflections calculated from the 
expanded deflection equations were only 0.017-0.529 percent off from those calculated from the 
exact deflection equations formulated for nonuniform cantilever beams. 

4. The further-developed displacement equations for nonuniform cantilever beams were used to 
formulate the displacement equations for nonuniform two-point supported beams and could be used 
for shape predictions of nonuniform structural panels. 

5. The proposed two-line strain-sensing system provided a powerful method for simultaneously 
monitoring the bending and torsional deformations of nonuniform cantilever beams such as wing 
boxes. The two-line strain-sensing system eliminated the need to install shear strain sensors to 
measure the surface distortions for monitoring cross-sectional rotations of the beam. 

6. A finite-element computer program was used to generate simulated strain data for input to the 
further-developed displacement equations. The displacement equations were then used to calculate 
deflections and cross-sectional twist angles, at discreet strain-sensing stations of the beam-like 
structure, for constructing the overall deformed shape of the whole structure. 

7. The high degree of accuracy of the further-developed displacement theory for nonunifornr beams 
was validated by finite-element analyses of various types of nonuniform beam structures (tapered 
tubular beams, such as cantilever and two-point supported beams; depth-tapered unswept and swept 
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wing boxes; width-tapered wing boxes; and double-tapered wing boxes) under combined bending 
and torsion. 


8. For tapered beams, the bending strain distribution along the axis of the beam was nonlinear, and the 
nonlinearity increased with increasing taper rate. The piecewise linear approximation of strain 
distributions was adequate to obtain highly accurate deflection curves for the whole range of beam 
taper rates. 

9. The Ko displacement theory, combined with the onboard strain-sensing system, form a powerful 
tool for in-flight deformed shape monitoring of unmanned aerospace vehicles (by ground-based 
pilots) to maintain safe flights. In addition, the real-time wing shape could be input to the aircraft 
control system for aeroelastic wing shape control. 
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FIGURES 



£(*) = £/_ 1 - (e w - £ ; ) X ^j = 1 ; < x < Xj 


Figure 1. Nonuniform cantilever beam instrumented with bending and distortion strain sensors. 



Figure 2. Deformed cantilever tubular beam showing slopes, dj , and deflections, >’/ , at strain-sensing 
stations, Xj(i = 0, 1,2,3, 
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1.0 



Strain-sensing station 


Figure 3. Relative magnitude of strain terms in deflection equation for a uniform cantilever beam. 



Figure 4. Linearly increasing downward shifting of deflection curve from cantilever beam position A to 
simple beam position B; zero shift at left end. 
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Figure 5. Tapered wing box instrumented with a two-line strain-sensing system. 



Figure 6. Tapered cantilever tubular beam instrumented with a two-line strain-sensing system. 
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Figure 7. Aircraft fuselage instrumented with a two-line strain-sensing system. 



Figure 8. Structural panel instrumented with multiple biaxial strain-sensing lines. 
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Deformed 


P= 100 lb 



/ =100.5 in. 
t = 0.02296 in. 
c 0 = 4 in. 
c n = 1, 2, 3, 4 in. 


Nodes 3673 

Four-node elements 3636 
Three-node elements 36 


090018 


(a) Tip vertical load of P= 100 lb; tip twisting moment of T= 0. 


P= 100 lb 



/ = 100.5 in. 
t = 0.02296 in. 
c 0 = 4 in. 
c n - 1 , 2, 3, 4 in. 


Nodes 3673 

Four-node elements 3636 
Three-node elements 36 


090019 


(b) Tip vertical load of P = 100 lb; tip twisting moment of T= 400 x c n in-lb. 


Figure 9. Finite-element model for a tapered cantilever tubular beam subjected to different tip loading 
conditions. 
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Figure 10. Bending strains calculated from SPAR element stress outputs for different tapered cantilever 
tubular beams; n = S;P= 100 lb. 



Figure 11. Comparison of deflections calculated from Ko displacement theory with those calculated from 
SPAR program for different tapered cantilever tubular beams; n = 8; P = 100 lb. 
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(a) Front strain-sensing line. 



(b) Rear strain-sensing line. 

Figure 12. SPAR-generated strains, { £;,£/ }, along front and rear strain-sensing lines on different tapered 
cantilever tubular beams; n = 8; P = 100 lb, T= 400 x c n in-lb. 


37 



(a) Front strain-sensing line. 



(b) Rear strain-sensing line. 

Figure 13. Comparison of deflections along strain-sensing lines calculated from Ko displacement theory 
with those calculated from SPAR program for different tapered cantilever tubular beams; n = 8; P = 100 lb, 
r=400xc„ in-lb. 


38 



Figure 14. Comparison of cross-sectional twist angles calculated from Ko displacement theory with those 
calculated from SPAR program for different tapered cantilever tubular beams; n = 8 ; P = 100 lb , 
T = 400 x c n in-lb . 
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/ =100.5 in. Nodes 3673 

t = 0.02296 in. Four-node elements 3636 

c 0 = 4 in. Three-node elements 36 

c„ = 1, 2, 3, 4 in. 

(a) Both ends simply supported. 


090027 



/ =100.5 in. 
t = 0.02296 in. 
c 0 = 4 in. 
c n = 1, 2, 3, 4 in. 


Nodes 3673 

Four-node elements 3636 
Three-node elements 36 


090028 


(b) Both ends fixed. 


Figure 15. Finite-element model for tapered tubular simple beam subjected to central downward vertical 
load of P = -200 lb. 
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Figure 16. Bending strains calculated from SPAR element stress outputs for different tapered tubular simple 
beams; n = 8; two ends simply supported; P = 100 lb. 


Strain-sensing station, i 



Figure 17. Comparison of deflections calculated from Ko displacement theory with those calculated from 
SPAR program for different tapered tubular simple beams; two ends simply supported; n = 8; P = 100 lb. 
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Strain-sensing station, i 

3 ’ 090031 

Figure 18. Bending strains calculated from SPAR element stress outputs for different tapered tubular simple 
beams; two ends fixed; n = 8; P = 100 lb. 


Strain-sensing station, i 



Figure 19. Comparison of deflections calculated from Ko displacement theory with those calculated from 
SPAR program for different tapered tubular simple beams; two ends fixed; n = 8; P = 100 lb. 
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-Depth-tapered unswept 
wing box (top view) 



/ = 100 in. 
w= 17 in. 
t = 0.02296 in. 


Nodes 3639 

Four-node elements 3632 
Three-node elements 10 


c 0 = 2 in. 

c n = 0.5, 1, 1.5, 2 in. 

(3 = 0 


090033 


(a) Unswept wing box. 



w=17in. Four-node elements 3632 

t = 0.02296 in. Three-node elements 10 


c 0 = 2 in. 

c n = 0.5, 1, 1.5, 2 in. 
(3 = 10 deg 


(b) Swept wing box. 

Figure 20. Finite-element models for depth-tapered wing boxes subjected to bending and torsion; 
P= 100 lb; T= 50 x win-lb. 
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(a) Front strain-sensing line. 



(b) Rear strain-sensing line. 

/ 

Figure 21. SPAR-generated strains, {£,,£, }, along front and rear strain-sensing lines on different depth- 
tapered unswept wing boxes; n = 8; P = 100 lb; T= 50 X w in-lb (clockwise). 
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(b) Rear strain-sensing line. 

Figure 22. Comparison of deflections, { )’i ,y \ }, along front and rear strain-sensing lines, calculated from 
Ko displacement theory with those calculated from SPAR program for different depth-tapered unswept wing 
boxes; n = 8; P = 100 lb; T = 50 x w in-lb (clockwise). 
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Figure 23. Comparison of cross-sectional twist angles calculated from Ko displacement theory with those 
calculated from SPAR program for different depth-tapered unswept wing boxes; n = 8; P = 100 lb; 

T= 50 x w in-lb (clockwise). 
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(a) Front strain-sensing line. 



(b) Rear strain-sensing line. 


Figure 24. SPAR-generated strains, {£/,£/}, along front and rear strain-sensing lines on different depth- 
tapered swept wing boxes; n = 8; P = 100 lb; T= 50 x w in-lb (clockwise). 
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1.0 

1.5 

2.0 


(a) Front strain-sensing line. 



Strain-sensing station, i 
(b) Rear strain-sensing line. 


1.0 

1.5 

2.0 


Figure 25. Comparison of deflections, { yj,y'i }, along front and rear strain-sensing lines calculated from Ko 
displacement theory with those calculated from SPAR program for different depth-tapered swept wing 
boxes; n = 8; P = 100 lb; T= 50 x w in-lb (clockwise). 
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Strain-sensing station, i 

3 090044 

Figure 26. Comparison of cross-sectional twist angles calculated from Ko displacement theory with those 
calculated from SPAR program for different depth-tapered swept wing boxes; n = 8; P = 100 lb; 
f=50xw in-lb (clockwise). 
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(b) Rear strain-sensing line. 


Figure 27. SPAR-generated strains, {£;,£/}, along front and rear strain-sensing lines on different depth- 
tapered swept wing boxes; n = 8; P = 100 lb; T = -(50 x w) in-lb (counterclockwise). 


50 


1.8 



Strain-sensing station, i 
(a) Front strain-sensing line. 


1.0 

1.5 

2.0 



Strain-sensing station, i 
(b) Rear strain-sensing line. 


1.0 

1.5 

2.0 


Figure 28. Comparison of deflections, { }’,■ , y\ }, along front and rear strain-sensing lines calculated from Ko 
displacement theory with those calculated from SPAR program for different depth-tapered swept wing 
boxes; n = 8; P = 100 lb; T = -(50 x w) in-lb (counterclockwise). 
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Strain-sensing station, i 

a 090049 

Figure 29. Comparison of cross-sectional twist angles calculated from Ko displacement theory with those 
calculated from SPAR program for different depth-tapered swept wing boxes; n = 8; P = 100 lb; 

T= -(50 x w) in-lb (counterclockwise). 


-Width-tapered 
wing box (top view) 


P= 100 lb 


/ =60 in. 
w 0 = 30 in. 



0.02296 in. 


c 0 = c n ~ 0.5 in. 
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Nodes 

Four-node elements 
Three-node elements 


3639 

3632 

10 


090050 


Figure 30. Finite-element model for a width-tapered unswept wing box subjected to bending and torsion; 
P = 100 lb; T= 50 x w n in-lb. 
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(a) Front strain-sensing line. 



Strain-sensing station, i 


(b) Rear strain-sensing line. 

Figure 31. SPAR-generated strains, {£/,£;}, along front and rear strain-sensing lines on different width- 
tapered unswept wing boxes; n = 8; P = 100 lb; T= 50 xw n in-lb (clockwise). 
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(a) Front strain- sensing line. 


10 

30 

20 



10 


20 

30 


(b) Rear strain-sensing line. 


Figure 32. Comparison of deflections, { yj,yj }, along front and rear strain-sensing lines calculated from Ko 
displacement theory with those calculated from SPAR program for different width-tapered unswept wing 
boxes; n = 8; P = 100 lb; T = 50 x w n in-lb (clockwise). 
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Figure 33. Comparison of cross-sectional twist angles calculated from Ko displacement theory with those 
calculated from SPAR program for different width-tapered unswept wing boxes; n = 8 ; P = 100 lb ; 
T= 50 xw n in- lb (clockwise). 
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Double-tapered 
wing box (top view) 


/ =100 in. 


w 0 = 34 in. 
w n = 17 in. 


0.02296 in. 



c 0 = 2 in. 
c n = 1 in. 


(3 =0 


Nodes 3639 

Four-node elements 3632 
Three-node elements 10 


Figure 34. Finite-element model for a double-tapered wing box subjected to bending and torsion; 
P = 100 lb, T= 50 x w n in- lb. 


x 10 -3 



Figure 35. SPAR-generated strains, { £/,£/ }, along front and rear strain-sensing lines on a double-tapered 
unswept wing box; n = 8; P = 100 lb; T= 50 x w n in-lb (clockwise). 
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Figure 36. Comparison of deflections, { y ; -,y z - }, along front and rear strain-sensing lines calculated from Ko 
displacement theory with those calculated from SPAR program, for a double-tapered unswept wing box; 
n = S;P = 100 lb; T= 50 xw n in-lb (clockwise). 



Strain-sensing station, i 


090059 


Figure 37. Comparison of cross-sectional twist angles calculated from Ko displacement theory with those 
calculated from SPAR program for a double-tapered unswept wing box; n = 8; P = 100 lb; T= 50 xw n in-lb 
(clockwise). 
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APPENDIX A 

DERIVATION OF DEFLECTION EQUATIONS FOR 
NONUNIFORM CANTILEVER BEAMS 

Slope equation (11) and deflection equation (13) for nonuniform cantilever beams are duplicated here as 
follows: 

Slope equation, 


tan 6 1 = Al 


£ i— 1 £ i . £ i-l c i £ i c i — 1 > c i 
c i—\~ c i (c ; -_i — Cj) c i— 1 


+ tan 6{_\ 


(Al) 


Deflection equation, 


J / = ( A/) 2 


£ i—l £ i 
2(c/-i _ c j ) 


£ i—\ c i £ i c i - 1 
( c i— 1 — c i ) 


Ci \og e 


c i 

c i - 1 


+ (c*-l 


D) 


+ V/-1 +4/tan0,_| 


(A2) 


Writing out equation (A2) for different indices, z, and making use of the indicial relationships expressed 
in equations (Al) and (A2), we obtain the following: 

For i = 1, 


VI = (Al) 2 


£ 0~ £ \ 
2(cq - q ) 


gpq ~ fwo 

(c 0 - q ) 3 


Cl 

q loge — + (c 0 -q) 
c o 


+ yo+ Al tan 6q 


(A2-a) 
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For i = 2, 


y 2 = (Al) 
= (Al) 


2 £ \ ~ £ 2 £ I c 2 ~ £ 2 C 1 


2 (q-Q) (q-c 2 ) 3 

2] £ 1 - £ 2 C|C 2 ~£ 2 c \ 


2(q -c 2 ) (q-c 2 ) 3 


c 2 lo g e — + (q-c 2 ) 

q 


c 2 l °S e — + (q - c 2) 


+ V| + Al tan 6\ 


q 


>’i 


+(Aiy 


£ 0 ” e l C()C| -qc 0 


2(co-q) (c 0 -q) 3 

Al tan 6\ 


c\ 

q lo §e — +(q>-q) 

c o 


+ jo + Al tan 0() 


+(Air 


£ o~ £ i t £ oq- £ i c o log q_ 

e c o 


c o-q (co-qr 


+ 4/ tan 6q 


= {Al) 


2 £ 1~ £ 2 £ 1 c 2 ~ e 2 q 


2(q -c 2 ) (q-c 2 ) 3 


c 2 


c 2 l°ge ~ + (q “ c 2 ) 


ft 


+ ■ 


gp - q £pq ~ £ i c o 


J 2(c 0 -q) (c 0 — q )~ 


Cl 

q ioge — + (co 

CO 


+ (AIY 


gp-q l gpci-gcp 1q q_ 

e c 0 


+ JO + 2 Al tan 6q 


co q (c 0 -q) 

After the terms are grouped, equation (A2-b) becomes 


> 2 = (Air 


£ 1~ £ 2 +3 _ £ 0 _£ 1 


2(q-c 2 ) 2(c 0 -q) 


(/U) 2 { £ l £ ' 2 - £ M 

Co 

c 2 log e — + (q - c 2 ) 

gpq-qcp 
' 0 

r q ij 

cilog e — + (c 0 -q) | 

. (q-c 2 ) 3 

L q J 

(cp - ci ) 

L c 0 jj 


+(Aiy 


£ ° q “ £lC °log e ^ 


(cp-q)“ 


cp 


+ )'o + 2 Al tan 9 q 



(A2-b) 


(A2-b') 
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For i = 3, 


V 3 = (A/) 2 - 

f £ 2- £ 3 

£ 2 C 3 - £ 3 C 2 

c 3 l°gg — 

L c 2 

[2(c 2 -c 3 ) 

(c 2 -c 3 ) 3 

= (AO 2 1 

£ 2 — £ 3 

£ 2 C 3 - e 3 c 2 

c 3 log e — 
c 2 

2(c 2 -c 3 ) 

(c 2 -c 3 ) 3 

f r 


r 



■ + }’2 + A/ tan 0 2 


V2 


(A/r 


£ 1 - £ 2 e l c 2~ e 2 c l 


2(q -c 2 ) 


(q -c' 2 ) L 


c 2 i°ge — Kq — c 2 ) 

q 


£ o~ £ i £ oq _£ i c o 


2 ( £ 0 -q) (c 0 -q) 3 L 


q !og e — + (c 0 -q) 
c 0 


+ (A If 


£ o- £ i A £ oq- £ i c o log q_ 
e c 0 


c o-q (c 0 -qr 


+ vo + 2 Al tan 6q 
A l tan 0 2 


+(a iy 


£ i ~ £ 2 ! £ i c 2 ~ £ 2 q lo „ c 2 ! £ o ~ £ i t £ oq ~ £ i c o lo „ q_ 
q ~ c 2 (q-c 2 ) 2 q c 0 ~q (c 0 -q) 2 e£>c o 


+ A / tan 9q 


= (A/) 


2 £ 2~ £ 3 £ 2 c 3 _£ 3 c 2 


2(c 2 -c 3 ) (c 2 -c 3 ) 3 L 

£ o - £ i £ oq - £ l c 0 


£ 3 lo ge — + (c 2 -c 3 ) 
c 2 


£ 1 _£ 2 £ l c 2 _£ 2 q 


2(q — c 2 ) (q-c 2 r L 


C'j 

^iogg— + (q -c 2 ) 
q 


2(co-q) (c 0 -q) 3 L 

2 


qioge — + (co-q) 

c 0 


+ (AZ) 


q^ 2 + £ F- 2 - £ 2 q loge £ 2 +2 ^I + 2 £ oq- £ i^Q i oge fL 
q ~ c 2 (q — c 2 ) 2 q c o-q (c 0 -q ) 2 c o 


+ yo + 3 Al tan 6q 


After the terms are grouped, equation (A2-c) becomes 


(A2-c) 


F3 = ( Al ) 
-( Al) 


£ 2~ £ 3 | 3 £ 1 ~ £ 2 | 5 £ 0“ £ 1 


2(c 2 -c 3 ) 2(q - c 2 ) 2(c 0 -q) 


£ 2 c 3 - e 3 c 2 

c 3 l«g£> — + (c 2 -c 3 ) 

£ i c 2 - £ 2 q 
1 ^ 

Co 

£ 2 lo gc— +(q-c 2 ) 

(c 2 -c 3 ) 3 

L c 2 J 

(q-c 2 r 

L q J 


+ 


£ oq - £ i £ o 


(co-q) L 


C i 

q !og e — + (co-q) 

£ 0 


+{Aiy 


£lC2 ” £2q io ge q + 2 £ oq- £ i £ o logg q 


(q -q) 


q (co-q) 


£ o 


+ JO + 3 Al tan 


(A2-c') 
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For i = 4, 


£ 3 - £ 4 

£3 c 4 -£ 4 c 3 

[2 (c 3 -c 4 ) 

( c 3 — c 4 ) 3 

f £ 3“ £ 4 

£ 3 £ 4 -£ 4 c 3 

}2(c 3 -c 4 ) 

(c 3 -c 4 ) 3 . 

£ 2- £ 3 
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2(c 2 -c 3 ) 
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, £ o - £ 1 

£ oq - £ i c o 
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(c 0 -q) 3 L 

£ l- £ 2 , £ I c 2 “ £ 2 c 1 
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1 
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£ 3 !og e — + (q-c 3 ) 


2(q -c 2 ) ( q — c' 2 )' 


c 2 l°ge *"(q -Q) + 


£ o ~ £ i £ oq ~ £ i c o 

2(q>-q) (c 0 -q) 3 


q log e — + (co-q) 
c 0 


2 £ 2~ £ 3 , £ 2 c 3 — e 3 c 2 , c 3 l0 £ l- £ 2 , o £ l c 2 “ £ 2q , „ c 2 , a £ 0 “ £ 1 , o £ 0q “ £ 1 £ 0 C 1 


-log e — + 2 


log e — +3- 


-log e — 


c 2 ~ c 3 (C2-C3) c 2 q~ c 2 (q — c 2 ) Cl c 0-q (c 0 -q) c o 

+ vo +4 A/ tan0n 

(A2-d) 


61 



After the terms are grouped, equation (A2-d) becomes 


)’4 = (MY 


e 3 _e 4 | g g 2 ~ e 3 | g £ 1 ~ g 2 | n £ 0 ~ £ 1 
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CO 
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( ./VZ - Q j 

Based on the indicial progression patterns in equations (A2-a)-(A2-d), the deflection, j 7 , can be written 
in a generalized form with two summations (with different summation limits) as 


y i = (M) 2 J J 
7=1 
1 — 1 

7=1 


e j~ 1 £ j £ j~l c j £ ] C ]~\ 


2(c j—\ ~ c ;) 


7-1 ( c j'-l — c j) 
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I', log, 
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7 


(cy-i-Cy) 


loge 


c 7-l 
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- + (c 7 _ 1 -c 7 ) 


C; 


7-1 


(A2-e) 


+ JO + (i)M tan 9q 


which is equation (15). 

Alternatively, based on the progression of indices in equations ( A2- a / )-(A2- d'), the generalized 
deflection equation can be rewritten in an alternative form as 


y i = (Aij 2 J j 
7=1 
i — 1 

\2 


[2(t-7) + l] £j 1 £j 


£ 7 '-l C 7 £ 7 C 7'-1 
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which is equation (16). 
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The general deflection equation (A2-f) is checked for accuracy using / = 4 as follows: 
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which agrees with equation (A2- d'). 


63 



APPENDIX B 

EXPANSION OF LOGARITHMIC TERMS 

This appendix discusses the mathematical processes of expanding the logarithmic terms (taken from 
Appendix E of reference 2). 

When the cantilever beam is slightly tapered (that is, c/ — > q_] ), the logarithmic terms in slope equation 
(14) and deflection equation (16) will approach zero (that is, log e (q/q_i)— >0). Therefore, the 
perturbation method must be used to expand the logarithmic term, log e (q/q_i) , in the neighborhood of 
(q/q_l ) ~ 1 to yield nonzero mathematical expressions, which are derived in the following sections. 

Expansion of log e (c ( . /Cj_, ) for Slope Equation 

When the beam is slightly tapered (that is, q— >q_i), the logarithmic term, log e (q/q_i) , can be 
expanded in the neighborhood of (q /q_] ) ~ 1 (that is, log e (q/q_i) = 0 ) (ref. 6) as follows: 
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2 

Carrying out the expansion and retaining up to the second-order terms in (q_| - q ) causes the 
logarithmic term, log e (q/q_i ), to take on the following form when q approaches q_j (q— >q_i): 
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Equation B2 is equation (17). 
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Expansion of [c,.Iog e (c ( ./c,._ 7 ) + (c jV/ -c ; )] for Deflection Equation 

When the beam is slightly tapered (that is, c/— >c/_i), the term, [c ; - log e (c//c,-_i) + (c ; -_i -c,-)] , in 
deflection equation (16) can be expanded in the neighborhood of (c ; /q_| ) ~ I (that is, log e (cj/cj_i) = 0). 
In light of the log e (cj/cj_i) expansion equation (Bl), the expansion, [c/ log e (c ; 7q_i) + (c,-_j -Cj)], can 

3 

be carried out up to the third-order terms in (c ; _| -c/) when c/ approaches c ; _j(Cj — > c;_] ) (ref. 6) as 
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Neglecting the fourth-order term in (c ; _| -c ; ) , we obtain the final expression for 
[ c i 1og e (fj / Ci_ i ) + ( i - c ; - )] when c ; approaches c/_j ( c, -> c;_! ) as 


C’ 

log e ' + (q_i Cj ) 

(c ; _i -q) 2 

r 3+ (c,-_i-c,)j 

L c i-l J 

6c, -1 

c i-l 


-^-^1 


c i-l 


(B4) 


which is equation (20). 
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APPENDIX C 

DERIVATIONS OF DEFLECTION EQUATION FOR 
SLIGHTLY NONUNIFORM CANTILEVER BEAMS 

Deflection equation (21) for the slightly nonuniform cantilever beam contains the slope term, 
Al tan0,_j , as follows: 
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In light of equation (18), the slope, tan0 M , can be written as 
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Writing out both equations (Cl) and (C2) for different indices (/ = 1, 2 , 3, 4, 5), we obtain 
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Based on equations (C3-a)-(C3-e), the deflection, y t , can be rewritten for different values of i as 
follows: 
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For i = 2, equation (C3-b) becomes 
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For i = 4, equation (C3-d) becomes 
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For i = 5 , equation (C 3 -e) becomes 
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In summary, 
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or. 
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Observing the indicial behavior in equations (C5-a)-(C5-e), we can establish a generalized displacement 
equation for any index, i (= 1, 2, 3, n) as 
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which is deflection equation (22). 
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Also, based on the indicial behavior in equations (C5-b')-(C5-e / ), the generalized deflection equation 
can be written in a more compact form with single summation as 
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(/= 1,2,3, 


which is deflection equation (23). 

The general deflection equation (C5-f) is checked for correctness by using i = 4. For i = 4, 
equation (C5-f) becomes 
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which agrees perfectly with equation (C5-d) and thus confirms the accuracy of general deflection equation 
(C5-f). 


Likewise, the alternative general deflection equation (C5-f / ) is also checked for accuracy by using i = 4. 
For / = 4, equation (C5- f 1 ') becomes 
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which agrees with equation (C5-d / ) and thus confirms the accuracy of the alternative general deflection 
equation (C5-d'). 
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APPENDIX D 

DERIVATIONS OF DEFLECTION EQUATIONS FOR 
UNIFORM CANTILEVER BEAMS 

When the nonuniformity of the beam diminishes (that is, c ( - = c / _ 1 —c, as in the uniform cantilever 
beam case), slope equations (18) and (19) degenerate into the following forms: 
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Furthermore, deflection equations (21) and (23) take on the following forms: 
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Writing out equation (D4) for indices, i = 1-5, we obtain 
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Based on these indicial progressions, deflection equations (D4-a)-(D4-e) can be written in generalized 
forms for the indices, i, as 
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which is identical to equation (28). 
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APPENDIX E 

RELATIVE MAGNITUDES OF STRAIN TERMS FOR 
UNIFORM CANTILEVER BEAMS 

For the study of relative magnitudes of the strain terms in uniform beam deflection equation (28), beam 
tip deflection equation (28f), duplicated as follows, is used: 
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The classical case of a uniform cantilever beam of a length, /, subjected to a tip vertical load, P , is 
considered in the analysis. For this case, the bending strain, £(x), is a linearly decreasing function of x, with 
a maximum at the fixed end and zero at the free end according to the expression 
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in which £ 0 is the maximum bending strain at the built-in end. 

The beam tip deflection, y n , for different values of n (= 1, 2, 4, 8) are provided in the remainder of this 
appendix. 

For n = 1, the strain at the beam tip is £± = 0 X £ 0 , and equation (El) yields 
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in which the strain, £ 0 , is expressed in terms of the moment, M( 0), using the strain-moment relationship 
(eq. (4)). 
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For n = 2, the strains at strain-sensing stations, i - 1-2, are £, =— £ 0 , £2 = — £() , and equation (El) 
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3 2 10 

For n = 4, the strains at strain-sensing stations, i = 1-4, are £ x = —£ 0 , £ 2 = — £ 0 , £, =— £ 0 ,£ 4 = —£ 0 , 
and equation (El) yields 4 4 4 4 
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